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SYNOPSIS 


The  Aluminum  Smelter  at  Kitimat  is  underlain  by  compressible  sedi¬ 
ments  to  a  maximum  depth  of  350  feet.  These  sediments  consist  mainly  of 
three  principal  strata,  gravel,  silt  and  clay.  The  buildings  are 
founded  on  a  thick,  dense  fill  placed  on  top  of  these  compressible 
soils.  Settlement  due  to  the  fill  and  building  loads  have  been  continu¬ 
ously  recorded  at  numerous  points  over  the  area  since  the  first  fill 
was  placed  in  1952.  At  the  time  of  writing  the  magnitude  of  settle¬ 
ment  had  exceeded  five  feet  in  some  axe  as  and  was  continuing  at  the 
rate  of  one  half  foot  per  year.  Some  structural  damage  to  the  long 
Potroom  buildings  appeared  to  be  the  result  of  differential  settlement 
and  accompanying  horizontal  movement  at  the  ground  surface. 

Theoretical  time  curves  which  match  the  observed  settlements 
were  obtained  by  computing  a  time  curve  for  each  of  the  three  principal 
layers.  From  these  time  curves  future  settlements  were  predicted. 
Coefficients  of  consolidation  and  compressive  indices  were  computed 
from  the  matchjuig  theoretical  curves.  These  coefficients  differed 
appreciably  from  those  obtained  by  conventional  laboratory  tests,  but 
agreed  with  values  which  could  be  expected  in  nature. 

Measurements  indicated  horizontal  movement  of  the  ground  surface 
toward  the  zone  of  maximum  settlement.  The  cause  of  this  movement  was 
shown  to  be  the  deflecting  beam  action  of  a  thick  surface  layer  of 

relatively  incompressible  soil  modified  somewhat  by  the  effect  of  con¬ 
solidation  in  the  direction  of  the  major  principal  stress  within  a 
deep  compressible  stratum. 
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THE  KITIMAT  SMELTER,  195k 


Potlines  1  and  2,  looking  Northwest 


CHAPTER  I 


INTRODUCTION 

The  Problem 

The  smeltersite  at  Kitimat  is  underlain  by  sedimentary  deposits 
of  a  complex  nature  including  sand,  gravel,  silt  and  marine  clay* 

These  sediments  are  normally  loaded  and  vary  in  thickness  to  a  maxi¬ 
mum  of  3?0  feet*  Fill,  varying  in  depth  to  a  maximum  of  30  feet,  was 
placed  over  an  area  of  approximately  1800  feet  by  1200  feet  to  bring 
the  working  site  to  a  suitable  elevation*  This  original  area  has 
since  been  considerably  enlarged*  Settlement  due  to  the  fill  load 
on  the  underlying  compressible  soil  has  been  continuously  recorded 
at  numerous  points  since  the  first  fill  was  placed  in  1952* 

On  June  28,  1957*  the  last  date  for  which  complete  records  are 
available,  the  total  maximum  settlement  had  amounted  to  feet* 

Of  this,  3  feet  had  taken  place  after  the  buildings  were  constructed* 
The  rate  of  settlement  was  proceeding  at  approximately  one  half  foot 
per  year  with  no  indication  of  any  change  in  the  immediate  future* 

Two  phases  of  the  settlement  problem  are  considered,  estimation  of 
rate  and  magnitude  of  future  settlements  from  past  records  and  the 
effect  of  settlement  on  the  buildings* 

The  plot  of  observed  settlement  versus  time  does  not  completely 
resemble  the  theoretical  time  curve  for  a  one-layered  consolidating 
system,  and  previous  attempts  to  explain  the  observed  settlement  on 
this  basis  have  not  been  entirely  successful*  This  study  attempts 
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to  explain  the  shape  of  the  observed  settlement-time  curve  by  con¬ 
sidering  the  effect  of  each  of  the  three  major  soil  types  contributing 
to  the  settlement*.  Using  the  assumptions  which  agree  with  the  past 
records,  the  future  settlements  are  predicted. 

Until  very  recently  structural  damage  due  to  settlement  was  not 
a  serious  problem  despite  the  fact  that  the  magnitude  of  vertical  move¬ 
ment  had  not  been  uniform  from  point  to  point  over  the  area.  During 
the  summer  of  1957  however,  many  building  columns  had  developed  a 
distinct  tilt.  Some  were  out  of  plumb  by  as  much  as  1§  inches  in 
10  feet  of  height  while  others  remained  vertical.  This  was  accompanied 
in  some  cases  by  cracking  of  the  concrete  pedestals  which  supported  the 
steel  columns.  A  portion  of  this  study  is  an  inquiry  into  the  role 
settlement  has  played  in  the  movement  of  these  columns. 

Geology 

The  geology  and  early  history  of  Kitimat  has  been  adequately 
described  in  a  recent  paper  by  R.  M.  Hardy  and  Charles  F.  Ripley  (9  )•* 
The  geological  conditions  are  typical  of  the  coastline  of  British 
Columbia  at  the  head  of  the  many  inlets  or  fjord  channels.  The 
smelt ersite  is  located  at  the  head  of  the  Kitimat  arm  of  Douglas 
Channel.  At  this  place  the  valley  is  approximately  2  miles  wide. 

The  mountain  slopes  rise  steeply  from  the  valley  sides.  It  is 
thought  that  this  valley  was  filled  with  sediments  during  the  glacial 
and  interglacial  periods  of  the  Pleistocene  age.  The  Kitimat  River  is 


Numbers  refer  to  bibliography  at  the  end  of  the  text. 
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entrenched  in  these  sediments  and  flows  southward  emptying  into  Douglas 
Channel  at  Kitimat*  A  high  ridge  approximately  3  miles  upstream  from 
the  smelter  seems  to  mark  the  terminus  of  one  glacial  period.  North¬ 
ward  from  this  ridge  the  land  is  from  300  to  500  feet  above  present 
sea  level*  Southward  from  the  ridge  the  land  slopes  gently  from 
approximate  elevation  50  to  the  sea*  The  colored  map.  Figure  1*,  is  a 
plan  of  the  Kitimat  Municipality,  and  shows  the  location  of  the 
smeltersite  with  respect  to  the  main  topographic  features* 

Beneath  the  portion  of  the  smeltersite  considered  in  this  report 
the  thickness  of  sediments  above  bedrock  varies  to  a  known  maximum  of 
350  feet*  The  geology  of  the  sediments  is  complex  because  they  include 
alluvial  fan  deposits  of  sand  and  gravel  laid  down  by  the  Kitimat 
River  running  south,  and  Moore  and  Anderson  Creeks  which  are  two 
streams  running  from  mountains  eastward  across  the  smeltersite* 

Deltaic  sand,  silt,  and  marine  clays  are  all  interlaid  in  an  irregular 
manner*  In  the  upper  20  feet  numerous  buried  logs,  stumps,  peat  and 
organic  silt  layers  have  been  found  indicating  considerable  forest 
growth  as  sediment  from  the  River  and  Creeks  was  being  deposited. 

Most  of  the  logs  were  remarkably  well  preserved*  Before  construction 
commenced  in  1951a  the  entire  area  was  covered  with  a  dense  forest 
growth*  An  upper  layer  of  from  about  3  to  15  feet  thick  was  composed 
chiefly  of  peat  and  organic  silt* 


* 


Figures  and  tables  are  grouped  together  following  the  written 
text* 
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From  a  grid  of  i^idely  spaced  deep  drill  holes  an  average  general¬ 
ized  soil  profile  was  obtained.  This  is  shown  in  Figure  2.  The  subsoil 
is  generally  composed  of  three  principal  layers  of  varying  thickness. 

The  upper  layer  is  chiefly  sand  and  gravel  with  silt  and  organic 
matter  irregularly  distributed  throughout.  The  middle  layer  is  chiefly 
nonplastic  silt  but  again  contains  stratifications  of  fine  sand  and 
organic  matter  throughout.  The  lower  layer  is  lean  clay  stratified  with 
silt.  The  water  table  was  very  near  the  ground  surface. 

Conventional  laboratory  tests  indicated  the  following  properties 
for  the  silt  and  clay  layers. ^ 

Silt  -  Specific  Gravity  2.72 

Natural  moisture  content  33  percent 
Generally  non  plastic 

Slope  of  e-log  P  curve  at  the  anticipated 
applied  load  O.lij. 

Clay  -  Specific  Gravity  2.72 

Natural  moisture  content  33  percent  generally 
close  to  or  above  the  liquid  limit 

Liquid  limit  28 

Plastic  limit  18 

Plasticity  index  10 

Slope  of  e-log  P  curve  at  the  anticipated 
applied  load  0.25 


*  Descriptions  of  conventional-  laboratory  soil  mechanics  tests 
can  be  found  in  Reference  15 • 
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Site  Development 

To  raise  the  working  surface  above  the  high  water  tides  and  to 
be  safe  from  flooding  by  mountain  freshets  from  Moore  and  ihderson 
Creeks,  a  fill  of  up  to  30  feet  in  depth  was  required  over  the  entire 
smeltersite#  Prior  to  placing  the  fill,  the  upper  layer  of  the 
organic  overburden  was  removed  from  the  site.  The  depth  of  stripping 
varied  from  approximately  3  to  15  feet.  Fill  was  placed  on  the  site 
immediately  following  the  stripping  operation.  This  material  which 
consisted  of  a  clean  sandy  gravel,  was  placed  in  thin  lifts  and 
compacted  with  a  cat-drawn  50- ton  pneumatic  tired  roller.  This  gave 
an  excellent  foundation  material  which  was  free  draining  and  dense. 

The  initial  area  filled,  and  the  one  of  concern  in  this  report, 
was  located  immediately  north  of  high  tide  mark  and  adjacent  to  the 
west  side  of  the  valley.  The  filling  operation  commenced  in  June 
1952,  and  was  completed  to  elevation  Ii5*5  above  sea  level  between 
September  and  November,  1952.  Figure  2  shows  the  net  applied  load 
on  this  area.  It  was  here  that  the  first  buildings.  Potlines  1  and  2, 
were  constructed  during  1953* 

South  of  this  area  fill  was  placed  on  the  tidal  flats  to  pro¬ 
vide  a  storage  area  for  the  raw  materials  unloaded  from  the  ships. 

This  fill  consisted  of  sandy  gravel  obtained  from  dredging  a  shipping 
channel  in  the  harbour  and  was  completed  to  elevation  28*5.  This  fill 
was  placed  between  May  and  November  1952. 
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North  of  Potlines  1  and  2,  other  areas  of  similar  size  have  been 
stripped  and  filled  in  a  similar  manner*  The  fill  at  the  site  of 
Potlines  3  to  5  was  placed  between  November  1953  and  October  195k* 

North  of  Moore  Creek  the  site  is  prepared  for  three  additional  pot¬ 
lines  which  when  completed  will  extend  a  total  distance  of  about  one 
half  a  mile  north  from  Potline  1* 

Settlement 

To  detemine  the  settlement  pattern  caused  by  the  fill  load, 

15  settlement  gauges  were  placed  in  the  fill  area  and  read  periodic¬ 
ally*  Each  gauge  consisted  of  a  3  by  3  foot  square  timber  base  to 
which  a  standard  2  inch  diameter  pipe  was  attached*  The  timber  base 
was  set  at  the  base  of  stripping  prior  to  the  placement  of  any  fill, 
and  the  2  inch  pipe  was  carried  through  the  fill  to  the  surface*  The 
settlement  was  measured  by  running  a  line  of  levels  to  the  gauge  from 
a  pemanent  bench  mark  on  bedrock  adjacent  to  the  site*  The  level  was 
a  precise  tilting  level  and  the  rod  was  read  to  the  nearest  0*005  feet* 
Short  circuits  and  frequent  checks  on  turning  points  resulted  in  very 
exact  work  considering  the  number  of  settlement  gauges,  the  distances, 
and  the  frequency  they  were  read*  It  is  felt  that  the  settlement 
readings  are  accurate  to  within  0*01  feet* 

The  bottom  one  foot  of  settlement  gauge  pipe  was  perforated,  and 
care  was  taken  to  insure  the  pipe  did  not  plug*  Since  the  completion 

of  the  fill  these  gauges  have  also  been  used  to  record  the  water  table. 
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This  is  done  by  lowering  a  bomb  carrying  two  bare  wire  ends  into  the 
gauge  pipe*  A  small  battery  feeds  current  into  the  wires  at  the  sur¬ 
face  while  an  ohmmeter  records  the  instant  the  bare  ends  make  contact 
with  the  water  in  the  pipe© 

After  the  buildings  were  constructed,  settlement  observations 
were  made  at  approximately  200  points  on  the  buildings#  The  observa¬ 
tion  points  are  specially  prepared  brass  plugs  set  in  the  concrete 
columns  and  are  placed  around  every  building  at  approximately  100 
foot  spacing#  Readings  are  made  periodically  using  the  same  method 
and  with  the  same  care  and  precision  as  for  reading  the  gauges# 

Figure  3  shows  a  plan  of  the  buildings  and  the  location  of  all  settle¬ 
ment  observation  points  in  this  area# 

The  magnitude  of  settlement  has  varied  considerably  over  the 
smelter  area  as  shown  by  Figure  h  which  is  a  plan  of  contours  of 
equal  total  settlement  to  June,  1957*  Plots  of  settlement  versus 
time  for  any  particular  settlement  gauge  show  important  differences  from 
the  general  shape  of  a  theoretical  time  curve.  This  is  especially  true 
for  the  early  stages  of  settlement#  The  observed  curve  plotted  to 
both  arithmetic  and  logarithmic  time  scale  shows  no  indication  of  any 
immediate  reduction  in  the  rate  of  settlement#  A  previous  attempt  to 
predict  the  amount  of  total  settlement  was  done  by  fitting  a  theoretical 
time  curve  to  the  observed  points  assuming  the  initial  irregularities 
were  similar  to  the  initial  settlement  from  a  conventional  laboratory 


consolidation  test* 
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To  attempt  to  gain  a  better  understanding  of  the  nature  of  con¬ 
solidation  within  the  compressible  soil  below  the  smeltersite,  theoreti¬ 
cal  time  curves  for  each  of  the  three  principal  layers  were  fitted  in 
such  a  way  that  the  sum  of  the  settlement  from  each  individual  curve 
was  equal  to  the  total  observed  settlement# 

The  Smelter  Buildings 

A  complete  unit  of  the  smelters ite  buildings  at  Kitimat  is  shown 
in  Figure  3#  The  first  to  be  constructed  were  Potlines  1  and  2,  and 
the  auxiliary  buildings  in  the  central  courtyard#  These  buildings 
covered  the  initial  fill  area#  Construction  commenced  in  the  fall  of 

1952  and  was  completed  by  August  1 95k  when  the  first  part  of  the 
smelter  went  into  production#  Fill  for  Potlines  3  to  5  was  placed  in 

1953  and  the  buildings  were  not  fully  completed  until  the  fall  of  1956# 
The  only  serious  structural  danage  which  had  developed  thus  far,  which 
could  be  attributed  to  settlement  occurred  in  Potlines  1  and  2#  These 
buildings  were  the  subject  of  -this  study# 

Typical  longitudinal  and  cross  sections  of  a  pot  room  are  shown  in 
Figures  $  and  6#  Each  Potroom  building  of  Potlines  1  and  2  are  essen¬ 
tially  structurally  similar#  They  are  long  relatively  flexible 
buildings  of  steel  frame  construction#  A  9  foot  concrete  column  rises 
from  a  spread  footing  founded  in  the  gravel  fill  and  supports  the  steel 
superstructure#  The  only  moving  equipment  within  the  potrooms  is  a 
20  ton  electric  crane  which  spans  the  width  of  the  building  and  runs 
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the  entire  length  on  a  steel  beam  22  feet  above  the  working  floor, 
this  floor  is  approximately  at  the  ground  surface.  Rigid  electrical 
connections,  the  operation  of  the  crane,  and  the  operation  of  each 
individual  pot  makes  differential  settlement  undesirable,  but  these 
limitations  have  not  as  yet  proved  to  be  as  serious  a  factor  as 
actual  damage  to  the  superstructure. 

Settlement  During  Construction 

The  first  footings  were  poured  almost  immediately  following 
completion  of  the  fill,  and  were  set  at  their  design  elevation.  It 
was  soon  realized  however  that  the  magnitude  of  settlement  would 
exceed  the  original  estimate,  and  would  require  longer  to  complete 
than  originally  planned.  Accordingly  all  footings  not  already  in 
place  were  set  at  an  elevation  so  adjusted  that  all  footings  would 
settle  to  a  level  position  by  January  1,  1955*  It  was  estimated  that 
settlement  occurring  after  this  date  would  not  cause  intolerable 
differential  movement. 

Between  pouring  of  the  concrete  footings  and  erecting  the 
steel  a  revised  estimate  indicated  still  more  expected  settlement.  To 
account  for  this,  and  to  still  provide  for  all  steel  columns  to  be 
level  by  January  1,  1955  9  the  concrete  columns  were  extended.  This 
double  field  adjustment  resulted  in  deviation  from  design  length  and 
elevations  by  up  to  one  foot  for  some  columns.  No  changes  were  made 

in  the  horizontal  directions.  A  complete  record  was  kept  of  all 
deviations  from  design  procedure.  Unfortunately  the  records  only 
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contain  recommended  changes*  Only  a  few  accurate  independent  checks 
were  made  in  the  field  to  determine  the  exact  setting  of  each  column, 
from  these  checks  it  is  assumed  that  they  were  set  precisely  as  recom¬ 
mended*  Other  work  in  the  field  has  shown  that  departures  from  the 
recommended  settings,  if  any,  were  small* 

Structural  Damage  Due  to  Settlement 

Following  construction  of  the  buildings,  settlement  continued  at 
almost  a  uniform  rate  until  the  summer  of  195 7 Figure  7  shows  con¬ 
tours  of  total  settlement  since  construction.  Settlement  previous  to 
this  time  would  have  no  bearing  on  any  structural  damage*  Until  1957 
only  minor  signs  of  damage  were  evident*  A  few  cracks  appeared  in  the 
concrete  block  walls,  and  a  few  diagonal  braces  were  slightly  buckled* 
This  in  no  way  contributed  to  structural  damage  or  unsightliness. 

By  the  summer  of  1957  however  a  serious  problem  had  arisen. 

Many  potroom  columns,  especially  on  the  west  end  had  developed  a  dis¬ 
tinctly  noticeable  lean.  The  column  tops  were  out  of  plumb  to  the  west 
by  as  much  as  lj  inches*  On  the  east  end  some  columns  leaned  eastward 
but  to  a  lesser  extent*  This  was  accompanied  by  cracking  of  the  con¬ 
crete  columns*  These  cracks  appeared  more  severe  at  columns  joined 

*  June  28,  1957*  the  last  date  for  which  a  complete  set  of  readings 
was  available  at  the  time  of  writing*  Later  readings,  where 
available,  were  included  in  this  report. 
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by  diagonal  bracing,  but  were  evident  at  other  columns  as  well.  Many 
of  the  diagonal  bracing  members  had  buckled  relieving  the  excess  stress, 
but  the  opposite  member  being  in  tension  and  having  no  way  to  relieve 
the  stress  had  in  some  cases  appeared  to  have  caused  the  bottom  of  the 
steel  column  to  shift  and  crack  its  concrete  support. 

Examination  of  the  expansion  joints  on  the  crane  rail  beam  showed 
that  all  except  a  very  few  on  the  east  end  appeared  closed  which  indica¬ 
ted  the  crane  rail  beam  was  in  compression. 

Some  of  tiie  more  severe  damage  was  repaired  immediately.  This 
was  accomplished  by  replacing  the  diagonal  bracing,  jacking  the  steel 
column,  and  chipping  out  and  replacing  the  damaged  concrete.  The  new 
diagonal  braces  were  fabricated  with  bolts  and  slotted  holes  so  that 
stress  could  be  periodically  relieved.  Less  severe  damage  was  left 
untouched  pending  the  results  of  an  investigation  to  determine  the 
primary  cause  of  damage  and  the  best  way  of  coping  with  the  problem. 

The  structural  damage  appeared  to  be  the  result  of  horizontal 
movements  caused  either  by  movement  of  the  ground  surface,  or  movement 
of  the  superstructure  due  to  some  external  force.  It  was  suggested  that 
perhaps  the  rapid  acceleration  and  stopping  of  the  electric  crane  dev¬ 
eloped  sufficient  horizontal  thrust  to  shift  the  network  of  steel 
columns.  Other  factors  indicated  that  the  source  of  movement  was 
within  the  soil  which  supported  the  footings. 
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Three  methods  of  approach  were  used  to  attempt  to  ascertain  and 
explain  any  horizontal  movements  within  the  soil*  These  were, 
theoretical  computation  of  horizontal  movement  due  to  consolidation, 
model  analysis,  and  comparison  to  the  problem  of  mining  subsidence* 


.  ,•  ’  •  ■  •  ■’  . 

:  •-  ;  ■  :  '  ■  -  "  ’  '  • 


CHAPTER  II 


MATCHING  THEORETICAL  TIME  CURVES  TO  OBSERVED  SETTLEMENTS 
Consolidation  Theory 

The  tern  consolidation,  as  used  in  soil  mechanics  explains  the 
settlement  which  occurs  within  a  mass  of  soil  subject  to  a  static  load* 
It  is  not  to  be  confused  with  mechanical  compaction  of  soil  under 
dynamic  loads,  or  to  solidification,  sometimes  called  consolidation 
by  geologists*  The  first  rigorous  mathematical  solution  to  this 
problem  was  published  by  Terzaghi  in  1923  and  is  found  summarized  in 
most  text  books  on  Soil  Mechanics  Stated  briefly,  the  Terzaghi  theory 
of  Consolidation  explains  the  settlement  which  occurs  due  to  an  applied 
load  on  a  layer  of  soil  as  being  a  result  of  the  soil  grains  squeezing 
closer  together*  This  reduces  the  void  spaces,  and  hence  the  volume  of 
the  soil,  resulting  in  settlement  at  the  surface* 

The  Terzaghi  theory  assumes  a  saturated  claylike  soil.  At  the 
instant  the  load  is  applied  it  is  entirely  carried  by  the  water  within 
the  voids  and  the  pressure  of  this  water  immediately  exceeds  hydro¬ 
static  pressure  by  the  amount  of  the  applied  load*  Under  this 
hydrostatic  excess  the  water  is  gradually  driven  from  the  void  spaces. 


See  for  example  Reference  2lp 
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the  soil  grains  move  closer  together  to  take  up  the  volume  once  occupied 
by  the  water,  the  load  is  gradually  transferred  from  the  water  to  the 
soil,  the  hydrostatic  excess  pressure  becomes  zero,  and  the  load 
carried  by  the  soil  grains  increases  by  an  amount  equal  to  the  load 
applied*  Thus  equilibrium  is  reached  under  this  new  loading  condition, 
and  the  process  of  consolidation  is  complete.  Any  additional  load 
change  would  cause  the  process  to  be  repeated. 

At  any  intermediate  time  the  degree  of  consolidation  can  be 
expressed  as  a  ratio  of  the  reduction  in  hydrostatic  excess  pressure  to 
the  initial  excess*  Since  reduction  in  hydrostatic  excess  pressure  is 
accompanied  by  volume  reduction  in  the  soil,  the  degree  of  consolidation 
can  also  be  expressed  as  the  ratio  of  settlement  at  any  time  to  the 
total  settlement  which  will  occur  under  the  particular  loading  condi¬ 
tions.  Because  settlements  are  generally  easier  to  observe  than  pore 
water  pressure  the  latter  is  the  more  common  way  of  expressing  degree 
of  consolidation.  The  mathematical  solution  to  the  problem  expresses 
the  degree  of  consolidation  U,  as  a  function  of  dimensionless  number  T, 
called  a  time  factor.  This  time  factor  depends  on  the  drainage  condi¬ 
tions  of  the  soil,  and  plots  of  time  factor  versus  degree  of  consolida¬ 
tion  are  available  for  several  drainage  conditions.  (21;  and  27) 

The  plot  of  time  factor  versus  the  average  degree  of  consolida¬ 
tion  within  a  layer  of  soil  is  called  a  theoretical  time  curve.  A 

normal  curve  begins  at  the  origin  and  ends  as  an  asymptote  to  the  100^ 
consolidation  line.  It  has  been  found  that  when  ‘the  time  factor  scale 
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is  changed  from  arithmetic  to  a  log  or  square  root  scale  a  large  portion 
of  the  curve  is  a  straight  line*  It  also  has  the  advantage  of  showing 
the  entire  curve  without  an  extremely  long  time  factor  scale.  The  use 
of  a  semi  log  plot  is  preferred  in  Western  Canada,  and  is  the  form 
adopted  for  this  study. 

The  actual  time  t,  required  for  a  consolidating  stratum,  of  soil 
to  reach  any  particular  degree  of  consolidation  is  determined  by  the 
equation : 

T  «  cv  t  *  (1) 


where  H  is  the  length  of  the  path  a  particle  of  water  must  travel  to 
the  nearest  drainage  face  from  the  farthest  point  from  this  face.  If 
the  stratum  is  drained  on  one  side  H  is  the  thickness  of  the  stratum. 

If  drained  on  both  sides  H  is  one  half  the  thickness.  The  term  cv  and 
H  are  essentially  constant  over  a  loading  increment  of  usual  magnitude. 

If  cv  and  H  are  known,  the  time  for  the  soil  to  reach  any  degree 
of  consolidation  under  an  applied  load  can  be  computed  by  equation  1. 

On  many  jobs  these  values  are  not  accurately  known.  If  records  of  the 
settlement  at  any  time  are  available  they  can  be  plotted,  and  by  suc¬ 
cessive  approximations  using  the  above  formula  a  theoretical  curve  can 


Subscripts  are  used  to  denote  a  particular  degree  of  consolidation. 

Thus  t^Q  and  T90  refer  to  the  time  at  which  9Q£  of  the  consolidation 
has  been  completed. 
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be  obtained  which  will  match  the  plot  of  observed  values.  This  process 
is  called  time  curve  fitting  or  matching.  Both  cv  and  H  are  assumed  to 
be  constant.  Also  by  this  method  a  theoretical  curve  can  be  fitted  to 
the  observed  values  before  the  process  of  consolidation  is  fully  com¬ 
plete.  Vihen  a  theoretical  curve  is  found  which  fits  the  observed 
settlement  it  can  then  be  extended  over  the  full  cycle  of  consolidation. 
Future  settlement  can  then  be  predicted  on  the  basis  that  it  follows 
the  fitted  theoretical  curve. 

The  first  problem  in  fitting  time  curves  is  to  choose  time 
factors  for  the  various  degrees  of  consolidation.  Taylor  (2lj.)  and 
others  have  published  time  factor  values  for  different  drainage  con¬ 
ditions  and  initial  excess  hydrostatic  pressures.  In/hen  these  values 
are  plotted  together  it  is  observed  that  all  curves  are  similar.  The 
choice  of  which  set  of  time  factors  to  use  is  difficult  because  the 
type  of  drainage  and  initial  excess  hydrostatic  pressure  are  usually 
not  accurately  known.  Taylor  gives  one  particular  theoretical  time 
curve  which  he  states  is  generally  accepted  as  an  adequate  representation 
of  typical  cases  in  nature.  It  is  for  one  dimensional  drainage  and 
unifom  initial  excess  hydrostatic  pressure.  Whether  or  not  a  case 
satisfies  all  the  assumptions  in  the  derivation  of  this  specific 
theoretical  curve,  especially  that  of  one  dimensional  drainage,  is 
open  to  question.  Unless  there  is  strong  evidence  to  support  the  use 
of  a  particular  theoretical  time  curve  then,  the  preferred  one  given  by 
Taylor  should  be  used.* 

*  Time  factor  values  for  this  case  are  given  in  the  Appendix. 
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The  Terzaghi  theory  of  consolidation  assumes  the  load  to  be 
applied  instantaneously,  and  time  factors  are  computed  on  this  basis* 

In  practise  this  simplification  is  seldom  realized*  The  load  is 
usually  applied  over  a  considerable  length  of  time  called  the  loading 
period*  A  method  of  correcting  theoretical  time  curves  for  this  loading 
period  has  been  advanced  by  Terzaghi  ( 2 5)*  This  method  assumes  that  the 
settlement,  or  degree  of  consolidation  at  the  end  of  the  period  of 
uniformly  increasing  load  would  be  the  same  as  if  the  total  load 
were  applied  suddenly  after  one  half  the  loading  period  had  elapsed* 

It  follows  that  if  P  is  'the  total  applied  load,  Pt  the  total 
load  applied  at  some  intermediate  time  within  the  loading  period  and 
U*  the  degree  of  consolidation  at  time  t/2  after  the  sudden  application 
of  Pj  then  the  total  degree  of  consolidation  at  time  t  is  equal  to 

u  >  u«  _£!  (2) 

p 

where  t  is  the  time  measured  from  the  start  of  the  loading  period,  and  U 
is  the  amount  of  final  settlement  due  to  consolidation  of  the  soil  under 
the  load  Pi*  At  any  time  following  completion  of  the  loading  period  the 
theoretical  settlement  can  be  computed  from  the  basic  equations  assuming 
the  total  load  was  applied  instantaneously  at  one  half  the  loading 
period* 
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Observed  Time  Curves 

A  test  frequently  run  on  soils  for  which  the  consolidation  charac¬ 
teristics  are  desired  is  called  the  Consolidation  Test**  This  test  is 
designed  to  duplicate  the  basic  assumptions  of  the  Terzaghi  theory  of 
consolidation*  In  the  laboratory  a  small  sample  of  soil  is  placed  in 
a  suitable  mold  where  it  is  free  to  drain  from  only  the  upper  and 
lower  surfaces*  A  load  is  applied  to  the  drainage  faces  and  the 
deformation  or  settlement  of  the  surface  of  the  soil  is  measured  at 
regular  time  intervals* 

Time  curves  from  a  laboratory  test  resemble  the  theoretical  time 
curve  in  all  respects  except  at  the  very  beginning,  and  at  the  end* 
Settlement  at  the  beginning  of  the  test  proceeds  at  a  much  faster  rate 
than  expected  from  a  theoretical  time  curve*  In  fact  it  appears  that 
some  settlement  occurs  almost  instantly*  This  initial  compression  is 
not  explained  by  the  Terzaghi  theory*  The  usual  explanation  for  the 
occurrence  of  initial  settlement  is  that  it  is  partly  due  to  the  com¬ 
pression  of  any  air  entrapped  in  the  soil.  (2it)  After  the  load  has 
remained  on  the  soil  for  a  considerable  length  of  time  the  laboratory 
settlement  curve  does  not  become  asymptotic  at  any  particular  amount 
of  compression,  but  continues  at  a  slow  rate  even  though  the  excess 


* 


See  Reference  l£  for  a  complete  description  of  this  and  other 
laboratory  soil  mechanics  tests* 
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hydrostatic  pressure  appears  to  be  zero.  This  phenomenon  is  also  not 
explained  by  Terzaghifs  theory  of  consolidation,  and  is  called 
secondary  compression.  It  is  probably  due  to  plastic  flow  or  gradual 
adjustment  of  the  soil  structure  under  the  new  load  (21;  and  26).  For 
most  inorganic  soils  both  initial  and  secondary  compression  are 
usually  small  compared  to  the  primary  compression  which  is  explained 
by  the  Terzaghi  theory.  No  proven  theory  exists  to  explain  the 
magnitude  or  rate  of  settlement  to  be  expected  from  initial  or 
secondary  compression.  The  question  of  whether  they  exist  in  the 
field  to  the  same  degree  as  in  the  laboratory  test  has  not  yet  been 
satisfactorily  answered.  Graphical  constructions  are  outlined  by 
Taylor  (21; )  to  determine  the  theoretical  zero  and  100  percent  con¬ 
solidation  from  a  laboratory  time  curve. 

Pressure  -  Void  ratio  Curves 

The  maximum  consolidation  pressure  to  which  a  soil  has  ever  been 
subjected  in  nature  is  called  its  preconsolidation  load.  When  this  load 
equals  the  present  overburden  pressure  the  soil  is  said  to  be  normally 
consolidated.  When  the  present  pressure  on  the  soil  is  less  than  it 
has  been  in  previous  geological  time  the  soil  is  said  to  be  precompres¬ 
sed  or  preconsolidated. 

When  a  consolidation  test  is  performed  on  a  sample  of  soil  it  is 
usually  loaded  in  small  increments  and  allowed  to  completely  con¬ 
solidate  under  each  increment  of  load  before  the  next  load  is  applied. 
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The  void  ratio  at  each  stage  of  the  test  can  be  calculated  from  the 
observed  compression  of  the  sample*  These  void  ratios  are  then  plot¬ 
ted  against  the  log  of  applied  pressure  and  the  resulting  curve  is 
called  an  e-log  P  curve,  or  a  pressure  void  ratio  curve.  For  pres¬ 
sures  in  excess  of  the  preconsolidation  pressure  the  e-log  P  curve  is 
usually  a  straight  line*  This  portion  of  the  curve  is  called  the  virgin 
branch.  The  virgin  branch  of  the  consolidation  curve  for  some  extra 
sensitive  soils  is  slightly  concave. 

The  slope  of  the  virgin  branch  of  the  e-log  P  curve  is  called 
the  compressive  index  (Gc)  of  the  soil.  Below  the  p  re  consolidation 
load,  or  where  the  e-log  P  curve  is  not  a  straight  line  settlement  is 
computed  on  the  basis  of  the  average  slope  of  the  curve  between  the 
pressures  of  interest. 

Knowing  the  compressive  index,  initial  void  ratio  (eQ)  and 
pressure  (p0)  of  a  soil,  the  void  ratio  (e)  after  any  pressure  change 
(  P),  can  be  computed  from  the  expression 

e  =  eQ  -  C0  Xogio  P0  P  (3) 

Po 

The  total  settlement  (S)  of  a  stratum  of  thickness  (H)  can  be  computed 
from 

S  s  e  "  eo  H*  (h) 


*  The  derivation  of  these  equations  can  be  found  in  any  text  book 
on  soil  mechanics,  for  example  see  Reference  27. 
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combining  these  two  equations  it  follows  that 

C0  =  S  (1  +  e0) 

H  iog10  P0  +  ^  P 

Po 


(5) 


Soil  Samples 

In  a  paper  entitled  Undisturbed  Clay  Samples  and  Undisturbed 
Clays  (26)  Terzaghi  compares  laboratory  consolidation  tests  to  the 
actual  field  behaviour  of  a  normally  loaded  clay  soil*  In  the  field, 
the  clay  has  been  deposited  by  sedimentation  at  an  estimated  rate  of 
O.Oij.  grams  per  square  centimeter  per  day*  This  rate  of  load  increase 
on  the  previously  deposited  soil  is  so  slow  that  certain  bond  stresses 
can  build  up  between  the  soil  grains.  These  bond  stresses  tend  to 
maintain  the  soil  in  a  state  similar  to  that  at  which  it  was  deposited. 
The  gradually  increasing  pressure  from  overlying  soil,  is  not  fast  or 
great  enough  to  destroy  these  bond  stresses  or  cause  the  amount  of 
compression  expected  from  laboratory  tests*  This  explains  why  the 
void  ratio  or  moisture  content  in  such  a  deposit  does  not  usually  show 
any  marked  decrease  with  increasing  depth. 

^Ihen  a  sample  of  soil  is  removed  from  the  ground,  the  total 
stress  on  it  is  decreased  to  atmospheric  pressure.  Capillary  forces 
however  tend  to  resist  any  volume  change*  $hen  tested  in  a  consolidation 
machine  the  load  is  reapplied  at  a  rate  of  up  to  several  tons  per  square 
centimeter  per  day*  This  greatly  exceeds  the  rate  of  loading  from 
material  deposition,  and  is  rapid  enough  to  destroy  the  bond  stress 
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between  the  individual  soil  grains*  With  these  bond  stresses  destroyed 
the  soil  sample  will  compress  to  a  smaller  volume  than  that  at  which  it 
existed  in  the  field  under  the  same  consolidation  pressure.  This 
explains  why  the  void  ratio  at  a  pressure  equal  to  the  preconsolidation 
load  determined  from  the  pressure  void  ratio  curve  is  smaller  than  the 
void  ratio  at  the  beginning  of  the  test*  Theoretically  capillary  forces 
should  hold  the  soil  from  expanding  when  removed  from  the  ground,  and 
if  no  destruction  of  the  internal  structure  of  the  soil  occurred  on 
reloading  there  should  be  no  change  in  void  ratio  for  loads  less  than 
the  preconsolidation  load*  Terzaghi,  as  just  stated,  attributes  this 
change  in  void  ratio  to  destruction  of  the  internal  structure  by 
rapid  loading*  Any  rough  handling,  or  change  of  orientation  of  the 
soil  grains  when  sampling  would  tend  to  further  aggravate  the  situation. 

The  rate  of  loading  in  the  field  from  man  made  sources  will 
approach  the  rate  of  loading  of  conventional  laboratory  tests*  It  is 
highly  probable  that  this  rate  of  field  loading  is  sufficient  to  at 
least  partially  destroy  the  internal  bond  stress  of  clay-like  soil.  For 
this  reason  compressive  index  values  from  laboratory  tests  may  not  be 
indicative  of  the  field  compressive  index  under  similar  loading  condi¬ 
tions.  A  method  of  determining  the  field  compressive  index  for  a 
normally  loaded  clay  soil  is  suggested  by  Terzaghi  and  Peck  (2? ) • 
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It  has  been  found  from  the  results  of  consolidation  tests  on 
soils  from  the  same  geological  area  that  all  of  the  virgin  compression 
branches,  if  extended,  intersect  at  approximately  the  zero  void  ratio 
axis.  Terzaghi  and  Peck  suggest  that  the  field  compression  curve  if 
extrapolated,  would  also  intersect  these  laboratory  curves  at  a  void 
ratio  of  zero*  This  then  is  one  point  on  the  field  virgin  compression 
branch*  The  existing  void  ratio  and  overburden  pressure  can  be  com¬ 
puted  from  a  few  relatively  undisturbed  samples  and  this  defines 
another  point  on  the  field  compression  curve*  Assuming  the  field 
virgin  branch  to  be  a  straight  line,  or  if  an  extra  sensitive  soil, 
slightly  concave,  the  field  virgin  compression  curve  can  be  drawn  and 
the  field  coup  res sive  index  determined* 

Load  Application  at  the  Smeltersite 

To  prepare  the  site  of  Potlines  1  and  2,  approximately  3  to  15 
feet  of  organic  soil  was  stripped  from  the  surface.  Stripping  was 
followed  almost  immediately  by  placing  of  the  fill*  This  fill  con¬ 
sisted  of  clean  sand  and  gravel  and  varied  in  final  thickness  to  a 
maximum  of  30  feet.  Shortly  after  completion  of  the  fill  to  its 
nominal  elevation  of  1*5*5  feet,  construction  of  the  buildings  commenced* 
A  few  footings  were  poured  late  in  1952,  but  the  major  construction 
occurred  in  1953*  Following  completion  of  the  buildings  additional 
fill  was  added  for  the  purpose  of  site  grading  and  landscaping. 

Figure  2  shows  a  plan  of  net  applied  load.  This  load  consists  of  the 
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weight  of  the  fill  plus  buildings  minus  the  weight  of  the  soil  stripped 
from  the  site.  The  weight  of  the  buildings  alone  amounted  to  only 
about  $00  pounds  per  square  foot* 

From  the  records  which  were  kept  at  the  time  the  following  is 
known  regarding  the  time  and  rate  of  application  of  these  loads. 

1.  The  date  each  settlement  gauge  was  installed  which  was  the 
approximate  date  of  placing  the  first  fill  in  the  vicinity 
of  the  gauge. 

2*  The  date  the  fill  was  completed  in  the  vicinity  of  each 
settlement  gauge. 

3*  The  date  the  site  grading  load  was  applied* 

it.  The  amount  of  the  site  grading  load  expressed  as  a  percent 
of  the  previously  applied  load. 

These  dates  and  magnitudes  are  shown  on  the  loading  diagram  of 
Figures  10  to  2f>. 

Bach  settlement  gauge  was  installed  prior  to  placing  any  fill  at 
the  particular  site  and  it  can  be  reasonably  assumed  that  the  total 
settlement  has  been  recorded  by  each  gauge.  The  actual  rate  of  load 
application  could  possibly  be  determined  from  daily  progress  reports* 
These  reports  were  not  readily  available ,  and  furthermore  the  loading 
period  was  so  short  compared  to  the  total  settlement  period  that  it 
was  felt  precise  determination  of  the  loading  curve  was  not  warranted* 
Examination  of  the  time  versus  settlement  plots  showed  very  little 
settlement  for  quite  some  time  following  installation  of  the  gauge. 
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This  was  followed  by  a  period  of  very  rapid  settlement  which  reduced 
rapidly  after  the  fill  was  completed.  The  time  settlement  curve  for 
each  gauge  was  examined  and  the  approximate  date  at  which  the  rapid 
settlement  commenced  was  taken  as  the  start  of  the  loading  period. 

This  date  is  shown  on  the  theoretical  time  curves.  Figures  10  to  25. 

The  load  was  assumed  to  increase  at  a  uniform  rate  from  this  date 
until  the  fill  operation  in  the  vicinity  of  the  gauge  was  completed. 

No  further  load  was  assumed  to  be  added  from  completion  of  the 
fill  until  the  completion  of  the  site  grading  load.  It  is  known  however, 
that  during  this  time  the  buildings  were  constructed,  but  their  weight 
was  small  in  comparison  to  the  weight  of  the  fill.  It  probably  took 
between  30  to  60  days  for  the  building  load  to  be  added.  The  site 
grading  load  was  assumed  to  be  instantly  applied.  In  reality  it 
probably  took  less  than  one  week. 

As  with  any  large  construction  operation  the  practise  is  seldom 
to  concentrate  on  one  specific  area  and  complete  it  before  starting  on 
the  next.  Most  likely  there  were  periods  of  comparative  inactivity 
within  each  assumed  loading  period,  these  xfould  probably  be  followed 
by  periods  of  very  intensive  loadings.  On  this  assumption,  and  based 
on  the  results  of  theoretical  time  curve  fitting  yet  to  be  described, 
a  probable  loading  diagram  at  each  gauge  was  shown  superimposed  over 
the  loading  diagram  just  described.  This  probable  loading  diagram  was 
not  essentially  different  from  the  first  assumed  loading,  and  was  not 
used  to  change  or  modify  any  of  the  computations  or  results. 
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Another  load  which  should  be  mentioned  is  increase  in  effective 
stress  resulting  from  a  lowering  of  the  ground  water  table.  Figures  8 
and  9  show  the  variation  of  the  ground  water  table  at  each  settlement 
gauge  since  1952.  Some  gauges  had  become  plugged  and  only  early  readings 
were  available.  Readings  between  195?  and  early  19?7  were  not  available 
at  the  time  of  writing.  The  water  table  at  each  gauge  has  lowered  con¬ 
siderably  since  the  fill  was  completed.  Most  of  this  lowering  however, 
appears  to  have  taken  place  within  the  first  two  years.  The  amount  of 
lowering  since  19?U  appears  to  have  been  generally  less  than  5  feet,  and 
may  be  within  the  range  of  seasonal  variation.  In  such  a  case  the 
effect  on  settlement  may  not  be  as  large  as  if  the  lowering  were 
permanent. 

Settlement  Observation 

Total  settlement  since  the  first  fill  was  placed  has  been  con¬ 
tinuously  recorded  by  means  of  1?  settlement  gauges  spread  over  the 
area.  One  piezometer  was  installed  adjacent  to  the  area  for  the 
purpose  of  observing  pore  water  pressures  within  the  consolidating 
soil.  This  installation  has  never  worked  properly  and  no  use  was 
made  of  any  of  the  readings.  The  1?  settlement  gauges  record  settle¬ 
ment  at  the  base  of  the  fill.  Settlement  observations  taken  simultane¬ 
ously  on  buildings  adjacent  to  these  gauges  have  shown  there  is  no 
measurable  settlement  within  the  fill  itself. 
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Settlement  Gauge  Time  Curves 

To  facilitate  the  computation  of  matching  theoretical  time 
curves  the  observed  settlement  was  plotted  to  a.  log  time  scale.  The 
settlement  at  each  gauge  was  plotted  on  a  separate  page  and  is  shown 
on  Figures  10  to  2£*  Figure  10  being  a  legend  for  these  settlement  time 
curves. 

Hatching  Theoretical  Time  Curves  to  Observed  Settlement 

It  has  been  shown  that  the  theoretical  tame  curve  for  a  single 
lawyer  of  soil  can  be  expressed  by  the  equation 
T  =  cv  t  (1) 


Assuming  cv  and  H  to  be  constant  over  the  loading  range  it  is  a  simple 
matter  to  compute  any  theoretical  time  curve.  By  trial  and  error  a 
theoretical  curve  may  be  found  which  will  match  settlement  observations 
providing  the  consolidation  is  proceeding  according  to  the  Terzaghi 
theory. 

The  sub  soil  profile  beneath  the  Kitimat  smelter  differs  con¬ 
siderably  from  the  ideal  assumptions  of  the  basic  consolidation  theory. 
The  most  important  difference  is  that  instead  of  a  one  layer  consoli¬ 
dating  system  there  are  three  principal  layers*  each  having  different 
consolidation  properties.  Besides  the  three  principal  layers  there 
are  countless  stratifications  and  differences  in  soil  properties 


within  each  stratum 
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Since  it  was  impossible  to  account  for  the  differences  within 
each  stratum  the  assumption  was  made  that  each  layer  was  in  itself 
homogenous,  isotropic,  and  in  all  respects  agreed  with  Terzaghi' s 
basic  assumptions  for  this  theory  of  consolidation.  This  of  course 
introduces  an  unknown  error,  but  it  is  the  only  known  way  of  dealing 
with  consolidation  problems,  and  is  the  method  nearly  always  used  when 
predicting  settlements  from  consolidation  tests# 

To  extend  the  Terzaghi  theory  for  a  one  layered  system  to  the 
three  consolidating  layers  at  Kitimat  a  theoretical  time  curve  was 
computed  for  each  layer:  The  first  problem  was  to  choose  a  proper 
set  of  time  factors.  It  was  decided  to  use  the  theoretical  time  curve 
published  by  Taylor  (2U)  for  linear  variation  of  initial  excess  hydro¬ 
static  stress.  The  basis  for  this  choice  was  that  all  theo retie al  curves 
are  quite  similar  and  as  Taylor  states,  "the  approximate  nature  of  the 
assumptions  in  the  Terzaghi  theory  generally  lead  to  the  accepted 
conclusion  that  this  curve  is  an  adequate  representation  of  typical 
cases  in  nature."  Time  factor  values  for  this  condition  are  found  in 
the  appendix. 

Having  chosen  the  theoretical  time  curve  it  was  a  straight 
f orward  matter  to  compute  a  specific  theoretical  curve  and  compare  it 
to  the  observed  plot.  This  was  done  by  means  of  equation  1. 
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To  begin,  a  certain  degree  of  consolidation  was  assumed  to  cor¬ 
respond  to  a  particular  time  and  settlement.  This  assumption  fixed  the 
entire  time  curve  which  was  then  easily  computed.  Equation  2  was  used 
to  correct  the  computed  theoretical  curve  for  the  settlement  during  the 
loading  period.  The  time  for  the  entire  curve  however,  was  reckoned 
from  the  beginning  rather  than  the  middle  of  the  loading  period.  The 
reasons  for  this  were  threefold.  First,  the  beginning  of  the  assumed 
loading  period  was  generally  quite  some  time  after  the  gauge  was  installed 
and  no  doubt  some  fill  had  been  placed  prior  to  the  beginning  of  the 
assumed  loading  period.  The  second  reason  was  that  the  loading  period 
was  always  less  -than  60  and  generally  less  than  30  days.  One  half  of 
this  time  would  not  make  any  appreciable  difference  in  the  log  time 
plot  after  the  completion  of  the  loading  period.  These  two  reasons 
coupled  with  idle  fact  that  the  load  was  assumed  to  be  applied  linearly 
when  in  reality  it  may  have  been  applied  at  some  other  rate  lead  to  the 
conclusion  that  it  mattered  little  whether  the  curve  was  computed  on 
the  basis  of  load  suddenly  applied  at  the  beginning,  or  middle  of  the 
loading  period. 

The  site  grading  load  was  accounted  for  by  assuming  that  the 
settlement  of  a  theoretical  time  curve  for  any  increment  of  load  would 
be  proportional  at  any  time  to  the  settlement  caused  by  any  other 
increment  of  applied  load.  Knowing  this  increase  in  load  compared  to 

the  initial  load,  and  the  time  this  increase  was  applied,  the 
theoretical  settlement  at  any  later  time  due  to  this  increased  load 
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was  computed  and  added  to  the  settlement  curve  for  the  initial  load. 
Since  the  computed  times  for  the  site  grading  load  did  not  generally 
coincide  with  the  computed  times  for  the  initial  curve,  the  addition 
of  settlement  was  done  by  plotting  the  theoretical  curve  for  one 
load,  and  then  graphically  adding  the  increased  settlement  due  to  the 
site  grading  load  at  the  proper  elapsed  time. 

By  this  method  a  theoretical  time  settlement  curve  for  each  of 
the  three  principal  layers  was  computed.  The  theoretical  total  settle¬ 
ment  of  the  entire  system  at  any  time  was  the  sum  of  the  settlements 
of  each  individual  layer.  Again  it  was  found  convenient  to  plot  the 
three  time  curves  together,  and  add  the  settlements  at  each  time 
graphically.  The  resulting  curve  joining  points  of  total  settlement 
for  different  times  was  the  theoretical  time  curve  for  the  entire 
system.  If  the  assumption  for  the  individual  curves  was  correct  then 
this  curve  matched  the  observed  settlement  time  plot.  If  the  curves 
did  not  agree  then  new  assumptions  were  made  for  each  individual  layer 
and  the  process  was  repeated.  Finally  by  trial  and  error  a  best  fit 
theoretical  curve  was  obtained. 

The  number  of  possible  combinations  is  countless.  It 
can  be  appreciated  that  this  process  was  tedious  and  time  con¬ 
suming.  However  with  practice  it  was  found  that  a  good  fit 
could  usually  be  obtained  with  comparative 
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ease.  It  required  approximately  three  quarters  of  an  hour  to  make  one 
complete  trial  and  seldom  was  a  good  fit  obtained  in  less  than  ten 
trials. 

The  best  fit  theoretical  time  curve  so  obtained  was  not  neces¬ 
sarily  a  unique  solution.  Other  assumptions  for  individual  layers,  when 
combined  could  also  give  a  time  curve  which  agreed  equally  well  with  the 
observed  plot.  For  this  reason  two  different  matching  theoretical  time 
curves  were  found  for  the  observed  settlement  plot  at  each  gauge.  To 
identify  the  two  fitted  curves  they  were  called  Series  A  and  Series  B. 
Series  A  was  the  fit  first  obtained  and  Series  B  was  the  second.  The 
results  are  shown  on  Figures  10  to  I4O  inclusive. 

After  Series  A  was  completed  it  was  noticed  that  the  assumed 
curve  for  the  gravel  layer  was  probably  almost  unique.  The  very  short 
period  of  consolidation,  and  the  difference  between  this  curve  and  the 
other  two  suggested  that  no  major  changes  should  be  attempted  for  this 
curve.  The  curves  for  the  silt  and  clay  layers  were  often  quite 
similar.  However,  at  many  of  the  gauges  one  of  the  curves  had  a  con¬ 
siderably  shorter  consolidation  period  than  the  other.  It  was  decided 
to  commence  the  fitting  for  Series  B  by  assuming  the  same  settlement 
for  each  layer,  but  interchanging  the  times.  The  gravel  layer  was 
generality  left  as  for  Series  A.  In  this  way  it  was  hoped  to  get  a  band 
of  possible  future  settlements  rather  than  just  one  theoretical  time 

settlement  curve  which  might  not  be  correct. 
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In  some  cases  this  interchange  of  times  had  no  effect  on  the 
final  total  theoretical  curve.  This  was  especially  true  in  cases  where 
the  assumed  times  for  each  layer  were  not  greatly  different.  At  other 
gauges  however*  the  mere  interchange  of  times  did  not  result  in  a  match¬ 
ing  curve.  In  these  cases  the  trial  and  error  procedure  was  again 
employed  to  furnish  a  curve  which  agreed  with  the  observed  one  and  yet 
was  composed  of  individual  curves  different  from  Series  A. 

The  results  of  these  two  separate  matchings  are  shown  compared 
on  Figures  lj.1  to  U8.  From  these  comparison  curves  it  can  be  seen  that 
the  only  significant  difference  between  the  two  methods  is  for  gauges 
11,  13  and  1U.  At  these  gauges  the  final  total  settlements  from  the 
two  fittings  differ  by  from  0.5  to  0.8  feet  or  by  approximately  10 
percent  of  the  total  settlement.  The  difference  in  theoretical  final 
total  settlement  at  all  other  gauges  was  less  than  0.20  feet. 

Identifying  Time  Curves  for  Each  Layer 

Having  completed  two  separate  fittings  for  each  gauge  there  still 
remained  the  question  of  which  single  curve  represented  the  settlement 
within  each  layer.  Gravel  soil  is  much  more  permeable  than  silt  or 
clay,  even  if  it  does  contain  considerable  organic  matter  and  other 
foreign  material.  It  should  therefore  consolidate  at  a  much  faster 
rate  than  either  the  silt  or  clay.  On  this  basis  the  single  curve 
showing  the  most  rapid  rate  of  consolidation  was  assumed  to  represent 
the  settlement  within  the  gravel  layer.  There  was  no  question  about 
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which  curve  to  choose  because  there  was  always  one  curve  which  reached 
approximately  90  percent  consolidation  about  the  time  the  loading  period 
was  completed*  Obviously  this  was  the  curve  for  the  gravel  layer* 

The  silt  and  clay  curves  were  more  difficult  to  determine  since 
they  were  often  not  greatly  different.  The  fact  that  a  good  match 
could  also  be  obtained  by  interchanging  the  times,  giving  Series  A  and 
Series  B  fits  further  complicated  'the  matter.  For  each  series  either 
individual  layer  curve  could  possibly  represent  the  consolidation  of 
either  the  silt  or  the  clay.  There  existed  four  possibilities: 

Series  A  with  the  middle  curve  representing  the  silt  and  lower  curve 
the  clay.  Series  A  with  the  middle  being  clay  and  the  lower  being  silt, 
and  Series  B  with  the  same  two  combinations.  The  choice  of  which  of  the 
four  represented  the  most  probable  combination  was  made  after  the 
coefficient  of  consolidation  and  compressive  index  were  computed  for 
each  possible  combination.  These  coefficients  were  based  on  results 
from  the  theoretical  time  curves  and  were  compared  to  values  from 
laboratory  tests. 

Coefficient  of  Consolidation 

From  the  results  of  several  laboratory  consolidation  tests  the 
coefficient  of  consolidation  was  computed  for  both  the  silt  and  clay 
soil  samples.  This  was  done  by  means  of  equation  1.  These  coefficients 
were  plotted  against  applied  load  as  shown  on  Figure  k9»  It  can  be 

seen  that  the  coefficient  of  consolidation  for  silt  is  approximately 
U  to  6  times  greater  than  for  clay.  The  range  of  values  for  silt  is 
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shown  to  be  from  approximately  1  to  3  ft*2  per  day  while  the  range  for 
clay  is  between  about  0*3  and  0*8  ft*2  per  day*  These  values  were  all 
computed  from  samples  which  were  loaded  and  drained  on  horizontal 
surfaces ,  that  is,  the  drainage  was  all  in  a  vertical  direction.  It 
was  appreciated  that  field  drainage  would  probably  tend  to  follow  any 
horizontal  lenses  of  high  permeability,  and  not  be  entirely  vertical* 
However  it  was  felt  that  these  laboratory  test  results  would  at  least 
give  some  indication  of  field  conditions,  and  relative  magnitudes  of 
cv  for  silt  compared  to  cv  for  clay* 

Equation  5  was  used  to  compute  the  field  coefficient  of  con¬ 
solidation.  Values  for  t  and  T  were  arbitrarily  chosen  for  the  condi¬ 
tion  of  90  percent  consolidation.  The  value  for  H  was  taken  as  the 
full  thickness  of  the  stratum.  This  assumed  drainage  occurred  from 
only  one  face  which  seemed  a  reasonable  assumption  considering  that 
the  lower  boundary  of  the  clay  is  bedrock  and  the  upper  boundary  of  the 
silt  is  gravel.  Certainly  the  silt  would  not  drain  into  the  clay,  but 
it  could  be  expected  that  the  clay  would  have  to  drain  into  the  silt, 
and  the  silt  into  the  gravel. 

Table  I  shows  the  data  and  results  of  the  computation  for 
coefficient  of  consolidation. 

Compressive  Index 

The  laboratory  compressive  index  or  average  slope  of  the 
e-log  P  curve  within  the  range  of  overburden  pressure  and  applied  load 
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averaged  O.II4.  for  the  silt  samples  and  0.25  for  the  clays.  The  shape 
of  the  e-log  P  curves  were  generally  quite  curved  indicating  that  the 
samples  had  suffered  some  disturbance.  Because  the  natural  moisture 
content  was  close  to  and  generally  above  the  liquid  limit  it  was  felt 
that  these  samples  probably  behaved  as  normally  loaded  clays  described 
in  Reference  26.  It  was  therefore  decided  to  attempt  to  obtain  a  field 
consolidation  curve  for  the  silt  and  clay  layers  by  the  method  outlined 
by  Terzaghi  and  Peck  (2?)«  As  previously  described  this  method  con¬ 
sisted  of  extending  the  laboratory  consolidation  curve  to  the  zero 
void  ratio  axis  and  joining  this  point  to  the  point  of  initial  void 
ratio  and  pressure. 

Accordingly,  a  number  of  e-log  P  curves  were  plotted  on  the  same 
sheet.  Trihen  these  curves  were  extrapolated  to  the  zero  void  ratio  axis 
it  was  found  that  they  all  intersected  very  close  to  the  same  point. 

By  joining  this  point  to  the  two  outside  limits  of  initial  condition 
a  band  of  probable  field  virgin  compression  curves  were  obtained. 

These  curves  are  shown  on  Figures  50  and  5l*  The  field  compressive 
indices  for  the  silt  ranged  from  0*25  to  O.I4J4.  and  for  the  clay  from 
0.3^  to  0.60. 

The  compressive  index  from  the  results  of  the  theoretical  time 
curves  were  computed  by  equation  5»  Before  this  equation  could  be 
applied  several  factors  had  to  be  evaluated.  The  settlement  (S)  was 

taken  directly  from  the  theoretical  time  curves.  The  thickness  of 
the  consolidating  stratum  (H)  was  taken  as  the  thickness  of  the  layer 
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as  shown  on  Figure  2.  To  obtain  the  initial  void  ratio  (e0)  an  examina¬ 
tion  was  made  from  the  available  laboratory  test  results.  Figure  £ 2 
shows  initial  void  ratios  computed  from  consolidation  tests  plotted 
against  depth  below  the  surface.  Terzaghi  (26)  explains  that  removing 
a  clay  like  sample  from  the  ground  does  not  appreciably  alter  the  void 
ratio.  It  was  therefore  assumed  that  the  initial  void  ratio  computed 
from  these  relatively  undisturbed  samples  was  probably  representative 
of  the  initial  void  ratio  existing  in  the  field. 

The  void  ratio  from  the  e-log  P  curve  at  the  overburden  pressure 
was  always  approximately  0.1  less  than  the  computed  initial  void  ratio. 
This  again  agrees  with  Terzaghi* s  hypothesis  for  normally  loaded  clays 


(26). 


Figure  £0  shows  that  the  initial  void  ratios  are  quite  scattered 
and  range  from  0.6  to  1.2.  The  majority  of  the  points  however  fall 
within  a  narrow  band  between  0.8  and  1.0  and  show  a  slight  decrease 
with  depth.  There  is  no  apparent  difference  between  silt  and  clay 
soils.  A  best  fit  straight  line  was  arbitrarily  drawn  through  the 
plotted  points  and  all  initial  void  ratios  for  future  calculations 
itfere  assumed  to  fall  on  this  line. 

The  unit  weight  of  the  soil  in  each  layer  was  also  required  in 
order  to  compute  the  overburden  pressure  (P0).  The  average  specific 
gravity  (Gs)  and  initial  void  ratio  of  both  the  silt  and  clay  was 

2.72  and  Q .90  respectively.  From  the  fundamental  equation 
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where  «  the  unit  weight  of  water,  the  buoyant  or  submerged  unit 
weight  (Y\))  s°il  was  computed  to  be  57  pounds  per  cubic  foot* 

This  unit  weight  was  used  in  computing  the  initial  stress  within  the 
soil. 

the  soil  within  the  gravel  layer  was  generally  a  coarse  dense 
gravel*  However  there  is  considerable  organic  matter  both  within  the 
voids  and  in  pockets  and  layers  throughout  the  entire  stratum*  The 
buoyant  unit  weight  of  the  gravel  was  arbitrarily  chosen  to  be  5 
pounds  per  cubic  foot  heavier  than  the  silt  or  the  clay  making  YD  £°r 
the  gravel  equal  to  62  pounds  per  cubic  foot.  Assuming  Gs  for  this 
material  to  be  equal  to  2*65  the  initial  void  ratio  was  calculated  by 
equation  6  to  be  0.66*  This  is  a  rather  low  unit  weight  and  high  void 
ratio  for  gravel  soil,  but  when  the  existence  of  organic  material  within 
the  stratum  was  considered  these  values  seemed  reasonable.  Using  these 
unit  weight  values  the  initial  stress  was  computed  at  the  mid  height  of 
each  principal  layer.  The  pertinent  data  and  results  are  shown  on 
Table  II. 

The  original  settlement  computation  done  by  the  engineering  con¬ 
sultants  assumed  an  initial  void  ratio  of  all  materials  to  be  equal  to 
1.0  at  a  pressure  of  1  ton  per  square  foot.  The  submerged  unit  weight 
was  assumed  to  be  60  pounds  per  cubic  foot  for  all  materials. 

Boussinesq  equations  and  Newmarks  influence  chart  (19) were  used 
to  compute  the  increase  in  vertical  stress  at  the  mid  height  of  each 
layer.  A  detailed  disciission  of  calculation  of  stresses  within  a  soil 
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due  to  a  load  applied  at  the  surface  is  given  in  Chapter  III.  Stiff  ice 
it  to  say  here  that  of  the  many  mathematical  solutions  to  this  type  of 
problem*  the  Boussinesq  solution  gives  the  largest  stresses*  and  his 
assumptions  probably  are  in  greater  disagreement  with  actual  soil  condi¬ 
tions  than  any  other  solution.  His  equations  were  used  however  because 
they  were  far  easier  to  apply  to  the  irregular  loading  conditions  than 
any  other*  and  they  are  commonly  used  by  others  in  the  field  of  soil 
mechanics.  The  results  of  this  stress  computation  are  shown  plotted 
against  depth  on  Figure  53  • 

Having  decided  upon  values  for  all  the  factors  the  compressive 
index  was  then  computed  from  equation  5*  The  data  and  results  of  this 
computation  for  all  the  combinations  are  shown  on  Tables  III  and  IV. 

Most  Probable  Curve 

"With  values  of  cv  and  Cccomputed  for  the  four  possible  combina¬ 
tions  using  theoretical  fitting  curves  Series  A  and  Series  B  the  most 
probable  theoretical  curve  could  then  be  dete mined.  This  was  done  by 
comparison  of  the  computed  values  with  results  from  laboratory  con¬ 
solidation  tests.  For  example*  the  possible  combinations  for 
Settlement  Gauge  1  were  as  follows: 
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SILT 

CLAY 

‘fcpQ  ~  l*i*oo 

t0Q  z  1*000 

Series  A,  middle  curve 

Assumed  to  be  silt, 

cv  =  5.55 

Cy  =  0.90 

Bottom  curve  assumed 
to  be  clay. 

s  =  1.60 

S  =  1.25 

C0  =  .221 

Cc  z  .766 

tpQ  ~  1*000 

tyQ  £  1*1*00 

Series  A,  middle  curve 

cv  z  6.10 

Cy  Z  0.82 

Assumed  to  be  clay- 
bottom  curve  assumed 

S  =  1.25 

S  z  1.60 

to  be  silt 

Cc  =  .172 

Cc  =  .98 

tj,0  z  l*ooo 

tpo  z  1*1*00 

Series  B,  middle  curve 

Cy  “  6.10 

Cy  =  0.82 

assumed  to  be  silt, 
bottom  curve  assumed 

S  r  1.60 

S  z  1.25 

to  be  clay. 

Cc  =  .221 

vO 

vQ 

C''- 

. 

1  1 

O 

o 

tpg  —  Wtoo 

tpo  —  1*000 

Series  B,  middle  curve 

Cy  =5.55 

cv  s  0.90 

assumed  to  be  clay, 
bottom  curve  assumed 

S  z  1.25 

S  s  1.60 

to  be  silt. 

cc  =  .172 

Cc  i  0.98 

The  units  were  t^o  in  days. 

S  in  feet  and  cv  in  f 

t. 2/ day.  ' 

The  first  of  these  four  possible  combinations  was 

assumed  to 

the  most  probable  field  condition. 

ho 

In  this  case  cv  and  Cc  for  the  silt  was  £.55  and  0.221  and  for 
the  clay  0.90  and  .766  respectively.  The  values  of  cv  for  silt  is 
somewhat  higher  than  computed  from  the  laboratory  tests.  The  silt 
compressive  index,  although  being  considerably  higher  than  that  obtained 
from  the  laboratory  tests,  was  only  slightly  lower  than  the  lower  limit 
of  the  range  of  field  compressive  inaicies  shown  on  Figure  $0.  The 
coefficient  of  consolidation  for  the  clay  agreed  with  the  results  of 
laboratory  tests.  The  compressive  index  was  however  slightly  higher 
than  the  highest  field  compressive  index  shown  on  Figure  $1.  The 
results  of  other  possible  combinations  showed  even  greater  disagreement 
with  laboratory  test  results  and  were  therefore  rejected. 

In  a  similar  manner  the  most  probable  set  of  time  curves  was 
chosen  for  each  settlement  gauge.  These  most  probable  values  for  cv 
and  Cc  are  shown  underlined  in  Tables  I,  III  and  IV.  They  are  also 
recorded  on  the  plot  of  the  most  probable  fitting  curve,  Figures  10  to 
25.  A  summary  of  the  results  of  theoretical  curve  matching  is  shown 
on  Tables  VI  and  VII. 

Significant  Figures 

All  of  the  computations  were  done  to  three  figure  accuracy.  The 
final  results  were  also  recorded  using  three  significant  figures.  It 
was  realized  however,  that  the  use  of  three  significant  figures  in 
computations  or  results  does  not  in  any  way  mean  that  these 
coefficients  were  known  to  this  accuracy.  In  consideration  of  the 
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uncertainty  of  the  many  assumptions  involved,  the  values  for  cv  and  Cc 
are  probably  not  even  accurate  to  one  significant  figure*  Three  figures 
are  required  to  insure  accuracy  in  the  mathematical  calculations,  and 
for  this  reason  only  are  three  significant  figures  shown  in  the  cal¬ 
culations  and  results*  These  physical  soil  properties  were  of 
secondary  importance  in  this  study,  being  used  mainly  to  identify  the 
most  probable  correct  time  curve*  Their  use  does  not  warrant  further 
discussion  on  the  accuracy  in  their  determination. 

Consolidation  Factors  for  the  Gravel  Stratum 

The  theoretical  time  curve  for  the  gravel  stratum  was  essentially 
the  same  for  both  Series  A  and  Series  B  fits.  Values  for  cv  and  Cc  were 
not  needed  to  assist  in  the  determination  of  which  curve  would  be  the 
most  probable  representation  of  field  conditions*  Computed  values  were 
used  only  to  obtain  an  indication  of  the  field  consolidation  character¬ 
istics.  It  was  appreciated  that  settlement  within  the  gravel  layer 
was  probably  mostly  a  result  of  decrease  in  volume  of  the  organic 
matter  contained  in  the  voids,  layers  and  pockets  of  this  stratum. 

Clean  dense  gravel  would  have  very  little  settlement. 

The  field  compressive  index  of  the  gravel  stratum  was  computed 
from  equation  5>.  The  factors  for  this  equation  were  obtained  in  a 
manner  similar  to  the  method  used  for  the  silt  and  clay.  The  settle¬ 
ment  data  was  taken  from  the  most  probable  fitting  theoretical  curve, 
and  no  alternative  values  were  computed.  These  compressive  indices 
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are  shown  computed  on  Table  V  and  are  summarized  on  Figures  11  to  25* 
and  Table  VII* 

No  attempt  was  made  to  compute  any  coefficients  of  consolidation 
for  the  gravel  layer*  It  was  felt  that  because  of  the  complex  nature  of 
these  surface  sediments  any  computed  cv  would  be  meaningless  if  not 
misleading* 

Initial  and  Secondary  Compression 

The  theoretical  time  curves  which  were  found  to  match  the 
observed  settlements  were  all  computed  on  the  basis  of  the  Terzaghi 
theory  of  consolidation.  As  previously  mentioned,  practical  experience 
has  found  that  the  total  settlement  of  a  soil  is  made  ip  not  only  of 
primary  settlement,  which  follows  the  theory  of  consolidation,  but 
also  of  initial  and  secondary  compressions*  These  latter  compressions 
do  not  follow  the  theory  of  consolidation  and  no  proven  theory  exists 
to  predict  the  rate  or  magnitude  of  settlement  from  either  cause.  No 
corrections  were  applied  to  the  matching  theoretical  curves  to  account 
for  any  initial  or  secondary  compression* 

Predicted  Future  Settlements 

The  purpose  of  obtaining  a  theoretical  time  curve  to  match  the 
observed  settlement  w as  to  explain  the  shape  of  the  observed  settlement 
curve,  and  predict  future  settlement  at  each  gauge*  The  matching 
Theoretical  curves.  Figures  11  to  1|0,  show  very  close  agreement 
with  the  observed  settlement*  Any  deviation  from  the 
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theoretical  curve  of  a  one  layered  system  consolidating  under  an  instant¬ 
aneously  applied  load  is  due  to  the  gradual  application  of  the  fill  load 
coupled  with  very  rapid  consolidation  within  the  gravel  layer.  This 
accounts  for  the  rapid  initial  settlement.  For  a  period  following  the 
completion  of  the  fill,  very  little  load  was  applied,  and  the  settle¬ 
ment  began  to  level  off  on  the  log  time  plot.  Increase  in  settlement 
on  the  log  time  plot  coincided  with  the  application  of  the  site  grading 
load  just  at  the  completion  of  building  construction.  I&thin  a  very  few 
days  following  the  application  of  this  load  settlement  within  the  gravel 
layer  had  reached  very  nearly  theoretical  100  percent  completion.  Small 
differences  of  up  to  approximately  0.3  feet  between  the  theoretical  and 
observed  curves  occur  before  this  time.  These  differences  are  probably 
due  to  deviations  from  the  assumed  actual  loading  conditions.  A  more 
probable  loading  diagram  is  shown  superimposed  on  the  assumed  diagram 
of  Figures  11  to  i;0. 

Following  placing  of  site  grading  fill,  the  seat  of  settlement 
was  almost  entirely  within  the  lower  silt  and  clay  layers.  This 
settlement  appears  to  be  following  basic  consolidation  theory  very 
nicely.  If  the  applied  load  remains  unchanged  future  settlements 
should  follow  the  fitted  theoretical  time  curve.  At  the  time  of  'die 
latest  set  of  readings,  September  1957,  the  total  settlement  varied 
from  60  to  85  percent  completed,  (Table  VI ).  The  settlement  writhin 
the  silt  and  clay  layers  was  between  50  and  70  percent  completed. 
According  to  these  fitted  time  curves,  this  remaining  settlement  will 
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take  place  at  a  gradually  decreasing  rate*  Most  of  this  settlement 
will  have  been  con^leted  within  the  next  ten  years,  and  theoretical 
settlement  beyond  I960  for  all  practical  purposes  will  be  zero.  Table 
VI  shows  a  summary  of  the  matching  theoretical  time  curves. 

If  however,  there  is  an  increase  in  applied  load,  the  magnitude 
and  rate  of  settlement  will  also  be  increased.  One  potential  cause  of 
increase  in  consolidation  load  is  the  increase  in  effective  stress  due 
to  lowering  of  the  ground  water  table.  Figures  7  and  8  show  the  latest 
ground  water  table  to  be  from  approximate  elevations  20  to  30*  These 
elevations  were  referenced  to  the  elevation  of  low  tide  water  which  is 
called  elevation  zero.  Since  high  tides  exceeding  20  feet  are  common 
in  the  ocean  immediately  adjacent  to  the  smeltersite  it  is  probable 
that  there  will  be  very  little  future  lowering  of  the  ground  water 
table. 

Another  factor  which  could  cause  the  settlement  to  exceed  that 
predicted  by  the  theoretical  time  curves  is  the  effect  of  secondary 
compression. 

How  much  secondary  compression  has  taken  place  in  the  gravel 
stratum  and  how  much  more  secondary  compression  will  occur  there  as 
well  as  in  the  silt  and  clay  layers  is  completely  unknown.  Further¬ 
more  there  is  no  known  method  of  estimating  this  settlement.  Judging 
from  laboratory  tests  on  inorganic  soils,  any  secondary  settlement  will 

probably  occur  at  a  very  slow  rate,  and  be  small  in  comparison  to  the 
primary  settlement. 
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Theoretical  Time  Curves  Compared 


Approximately  two  years  before  this  writing  a  different  theoreti¬ 
cal  time  curve  was  fitted  to  the  observed  settlement  for  gauges  1  to  15* 
This  fitting  was  done  by  the  conventional  method  of  assuming  only  one 
layer  of  consolidating  soil.  The  irregularities  in  the  early  part  of 
the  observed  time  settlement  pattern  were  treated  as  if  they  were  initial 
compression  in  an  ordinary  consolidation  test.  By  this  method  a  theoreti¬ 
cal  time  curve  was  obtained  which  matched  the  later  observed  settlements 
as  well  as  the  theoretical  curves  obtained  from  consideration  of  a 
three  layered  system.  Figures  i;l  to  i|8  show  a  comparison  of  the  time 
curves  obtained  from  this  one  layered  fitting  and  the  two  curves  from 
the  three  layered  consideration.  The  agreement  is  remarkable. 

Generally  the  theoretical  total  settlement  is  within  a  few  tenths  of  a 
foot.  The  only  exceptions  are  gauges  5  and  lh  which  disagree  in  final 
total  settlement  by  approximately  one  foot. 

A  study  of  the  coirparisons  between  theoretical  time  curves 
obtained  from  these  different  assumptions  leads  to  the  conclusion  that 
a  satisfactory  fit  can  be  obtained  from  the  simple  assumption  of  a  one 
layered  consolidating  system.  Certainly  it  is  much  faster  and  easier 
to  match  one  curve  to  observed  settlements  than  to  match  the  combined 
sum  of  three  separate  theoretical  curves.  There  is  little  to  choose 
between  the  final  end  result.  The  one  layered  curve,  however,  reveals 
nothing  of  the  behaviour  of  the  individual  soil  layers.  If  this 
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infoimation  is  required  then  the  more  tedious  three  layered  procedure 
must  be  employe d. 

Coefficient  of  Consolidation,  and  Compressive  Index 

A  by-product  of  the  theoretical  time  curve  fitting,  was  the 
determination  of  values  for  field  coefficients  of  consolidation  and 
compressive  indices  at  each  settlement  gauge.  The  main  assumptions 
and  methods  of  computation  have  already  been  discussed.  It  should  how¬ 
ever  be  repeated  that  the  value  for  these  consolidation  factors  are  no 
more  accurate  than  the  matching  theoretical  time  curves,  and  the 
assumed  soil  properties.  Although  they  are  recorded  to  three  figure 
accuracy  for  convenience  in  calculating  and  checking  the  mathematics 
they  are  probably  not  accurate  to  more  than  one  significant  figure. 

A  summary  of  cv  and  Cc  values  is  shown  on  Table  VII. 

These  consolidation  factors  were  found  to  vary  widely  from  one 
gauge  to  the  next.  Although  the  values  differed  considerably  from 
those  obtained  by  conventional  laboratory  tests  these  differences  can 
be  explained  by  consideration  of  the  physical  conditions  of  the  test 
compared  to  that  in  nature. 

The  coefficient  of  consolidation  for  both  the  silt  and  clay 
was  often  3  and  sometimes  as  much  as  10  times  that  obtained  by  labora¬ 
tory  tests.  The  laboratory  test  however,  pemits  drainage  in  only  the 
vertical  direction  while  in  the  field  drainage  will  occur  along  the 
path  of  greatest  permeability.  A  well  known  fact,  which  is  supported 
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by  sound  theory  (25)  is  that  for  a  stratified  soil  the  coefficient  of 
permeability  is  always  greater  parallel  to  the  stratifications  than 
perpendicular  to  them.  The  coefficient  of  consolidation,  being  pro¬ 
portional  to  the  coefficient  of  permeability  would  also  be  greatest 
in  a  direction  parallel  to  the  stratifications.  It  is  expected  that 
since  the  silts  and  clays  below  the  Kitimat  smelter  were  deposited  as 
marine  sediments,  they  would  contain  horizontal  stratifications.  This 
assumption  is  verified  in  part  by  the  drill  hole  logs  which  record  layers 
or  pockets  of  fine  sand  and  silt  within  these  strata.  It  is  expected 
that  drainage  of  the  silt  and  clay  layers  would  at  least  partly  follow 
these  horizontal  permeable  layers  and  thus  effectively  increase  the 
overall  coefficient  of  consolidation  from  that  of  vertical  drainage 
only. 

In  a  study  of  the  field  consolidation  characteristics  of  a 
uniform  deposit  of  Boston  blue  clay  which  to  visual  inspection  showed 
only  a  few  fine  sand  strata,  Gould  (8)  detennined  the  horizontal 
permeability  to  be  from  3h  to  50  times  the  vertical.  He  also  found 
for  this  material  that  the  time  rate  of  consolidation  from  three 
dimensional  drainage  was  about  3  times  that  expected  from  basic  theory 
of  vertical  drainage.  The  coefficients  of  consolidation  for  the  Kitimat 
sediments  computed  from  the  matching  time  curves  seem  quite  reasonable 
when  considered  in  the  light  of  probable  field  conditions. 
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Regarding  the  compressive  indices  it  was  noted  that  they  were 
considerably  larger  than  expected  from  laboratory  test.  Of  course  no 
tests  were  run  on  the  soil  within  the  gravel  layer*  When  gravel  is 
washed  into  place  by  a  river  or  stream  the  deposit  is  usually  quite 
dense  and  little  settlement  is  expected  from  a  dense  gravel  deposit* 

The  presence  of  compressible  material  within  the  gravel  contributes  to 
increasing  the  settlement  within  the  layer*  The  compressive  indices 
computed  from  the  matching  theoretical  time  curves  ranged  from  0*01 
to  0*12  with  a  general  average  of  about  0.06*  Since  it  is  almost 
impossible  to  make  laboratory  consolidation  tests  on  such  a  heterogeneous 
gravel  deposit,  about  the  only  way  of  determining  the  consolidation 
characteristics  is  by  field  observations  such  as  have  been  done  here* 

The  laboratory  compressive  indices  of  the  silts  and  clays  were 
taken  as  the  average  slope  of  the  e-log  P  curve  within  the  range  of 
field  loading*  These  values  were  considerably  lower  than  those  com¬ 
puted  as  shown  on  Figures  $ 0  and  5>1  which  agreed  quite  well  with  the 
time  curve  computations*  Some  of  the  indices  for  silt  were  slightly  low, 
and  some  values  for  clay  were  slightly  high*  The  general  agreement  with 
the  values  from  Figures  5>0  and  5l,  and  the  disagreement  with  the  values 
from  conventional  laboratory  tests  appeared  quite  significant*  It  would 
seem  that  for  a  normally  loaded  soil  of  this  type,  the  more  realistic 
procedure  for  estimating  settlement  is  that  based  on  a  field  compres¬ 
sive  index  determined  by  the  method  outlined  on  Figures  $0  and  5 1  rather 
than  from  the  average  slope  of  the  e-log  P  curve  over  the  range  of 
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applied  load.  This  conclusion  is  in  agreement  with  Terzaghi»s  hypothesis 
concerning  the  effect  of  sampling  a  normally  consolidated  clay  soil.  (26) 

Conclusions 

In  conclusion  it  can  be  said  that  a  theoretical  time  curve  was 
found  to  match  the  observed  settlement  at  each  gauge.  This  curve  was 
the  summation  of  settlements  from  time  curves  for  each  of  three 
principal  layers.  From  time  curves  of  each  layer  values  of  cv  and  Cc 
were  computed  for  each  soil  type.  These  consolidation  factors  were 
somewhat  different  from  those  expected  from  basic  theory  and  conven¬ 
tional  laboratory  tests,  but  agreed  quite  well  with  that  which  should 
be  expected  in  nature.  A  method  of  computing  the  field  compressive 
index  from  several  consolidation  tests  was  verified  for  the  Kitimat 
silts  and  clays. 

The  rate  and  magnitude  of  settlement  was  such  that  between  60 
and  80  percent  of  the  final  total  settlement  had  been  completed  by 
September  1957*  Settlement  within  the  gravel  layer  was  completed  very 
soon  after  the  buildings  were  constructed.  If*  the  loading  remains 
unchanged,  all  settlement  after  October  1957  will  be  within  the  lower 
silt  and  clay  strata  and  will  continue  at  an  ever  decreasing  rate 
until  after  1970. 
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CHAPTER  III 


CONSOLIDATION  AND  HORIZONTAL  DISPLACEMENT 


Theory 

Any  element  of  soil  below  a  surface  applied  load  is  acted  on  by 
stresses  in  all  directions.  In  the  Cartesian  system  of  coordinates 
six  independent  stresses  are  both  necessary  and  sufficient  to  completely 
define  the  stress  on  any  element*  (28)  In  dealing  with  settlement 
problems  the  vertical  component  of  stress  is  usually  all  that  is  of 
interest  because  the  settlement  or  movement  in  a  vertical  direction  is 
generally  the  only  concern*  Consolidation,  being  a  direct  function  of 
applied  excess  stress  will  be  a  maximum  in  the  direction  of  the  maximum 
stress  on  the  element  which  is  called  the  major  principal  stress.  Only 
directly  under  the  center  of  a  uniformly  distributed  load,  or  point 
load  can  the  direction  of  the  major  principal  stress  be  assumed  vertical 
without  detailed  analysis.  Since  the  applied  load  was  large  and 
irregular  it  was  decided  to  determine  the  direction  of  the  major  princi¬ 
pal  stress  at  various  points  within  the  consolidating  soil  beneath 
Potlines  1  and  2.  Knowing  this  stress  the  theoretical  horizontal 
components  of  movement  could  then  be  calculated  on  the  basis  that  the 
maximum  movement  was  in  the  direction  of  the  major  principal  stress. 

This  information  could  then  be  compared  to  the  actual  observed  move¬ 
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Stress  Computations 

A  search  of  the  pertinent  literature  showed  that  there  were 
three  distinctly  separate  methods  of  computing  the  stresses  in  an 
elastic  system  below  a  loaded  area.  These  methods  are  named  after  the 
men  who  first  advanced  the  theory^  Boussinesq,  Westergaard,  and 
Bumister. 

A  brief  review  of  the  assumptions  and  theories  of  each  of 
these  solutions  follows. 

Boussinesq  Equations 

The  first  adequate  solution  to  the  problem  of  stresses  within  an 
elastic  medium  acted  on  by  a  force  at  its  surface  was  published  by 
Boussinesq  in  1885*  His  equations  first  appeared  in  Gartesian 
coordinates  (30)  but  are  quoted  in  most  textbooks  on  soil  mechanics 
using  spherical  coordinates  (2l|).  Boussinesq  assumed  a  semi-infinite 
elastic  medium  which  was  homogeneous  and  isotropic.  His  equations 
define  all  components  of  stress  on  any  element  beneath  a  point  load 
acting  normal  to  the  surface.  The  derivation  of  these  equations  is 
found  in  several  places  (18,  28).  Later  authors  have  extended  his  work 
by  integrating  the  equations  for  point  loading  over  an  area,  and  so 
found  stress  equations  due  to  various  loading  conditions.  Other 
authors  have  modified  his  equations  to  account  for  departures  from  the 
original  assumptions  of  a  homogeneous,  isotropic,  semi-infinite 
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elastic  solid*  The  list  of  such  special  considerations  was  long  but  a 
few  however  were  worth  special  mention* 

In  192k,  Carothers  published  formulae  for  line  and  strip  loads, 
for  both  a  semi-infinite  mass,  and  for  the  case  of  a  smooth  and  rough 
rigid  boundary  at  some  finite  depth  below  the  surface.  The  solution 
of  the  case  of  a  rough  boundary  at  a  finite  depth  contained  an  error 
which  was  objected  to  at  the  time  in  an  oral  discussion  by  Timmoshenko, 
but  was  not  generally  known  until  years  later.* 

In  193U*  Jurgenson  (12)  expanded  Carothers1  equations  and  pub¬ 
lished  tables  and  charts  to  facilitate  the  computation.  These  charts 
and  tables  are  for  stresses  caused  by  line,  strip  or  symmetrical  area 
loads.  They  give  equations  in  one  plane  only.  Carothers'  error  is 
inherent  in  Jurgenson' s  solution  for  the  one  particular  case.* 

Biot  (1)  in  1935  derived  formulae  for  stress  at  the  rigid 
boundary  due  to  a  point  load  acting  on  the  surface  and  showed  that  the 
stress  so  computed  is  different  from  that  obtained  by  the  Boussinesq 
solution.  Cummings  (5)  expanded  Biot's  formula  to  apply  to  the  case  of 
a  circular  area  and  showed  that  the  stress  so  obtained  at  the  rigid 
boundary  could  be  approximated  by  using  a  depth  of  0.75  tines  'the 
actual  depth  in  the  Boussinesq  equations.  He  showed  that  the  stress 
on  the  rigid  boundary  is  greater  than  the  stress  obtained  from 

The  reference  for  this  is  Pickett  (21)  Page  35*  He  refers  to 
Carothers,  S.D.,  Elastic  Equivalence  of  Statically  Equipollent 
loads,  Proceedings  of  the  International  Mathematical  Congress, 

Vol.  2.  Toronto  University  Press. 
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Boussinesq^  solution  but  suggests  that  by  St.  Venants  Principal  the 
stress  concentration  would  be  relieved  a  short  distance  from  the  rigid 
boundary  and  be  the  same  as  computed  by  Boussinesq. 

In  1938,  Pickett  (21)  published  the  necessary  equations  for 
determination  of  the  stresses,  vertical,  horizontal,  and  shearing,  at 
any  point  in  a  homogeneous,  isotropic  elastic  solid  bounded  on  the 
lower  side  by  a  rough  rigid  boundary  and  carrying  a  known  normal  load 
on  the  surface.  This  was  complete  with  a  description  of  how  to  solve 
the  equations  for  special  cases.  The  solution  and  application  to  a 
specific  case  however  is  quite  involved.  He  states  that  for  vertical 
stress  the  existing  Newmark  Charts  (19)  are  satisfactory  for  Poisson1  s 
ratio  between  0.25  and  0.5* 

These  and  other  papers  on  variations  of  the  Boussinesq  equations 
were  generally  unsuitable  for  this  study  because: 

1.  The  equations  were  left  in  a  purely  mathematical  form,  not 
sufficiently  simplified  to  be  readily  usable. 

2.  Miere  simplified  formulae,  tables  or  charts  were  published 
to  facilitate  the  computations  they  dealt  only  with  stress 
components  in  one  plane,  or  with  very  simple  loaded  areas. 

3.  The  charts  and  tables  available  for  simplified  computations 
applied  only  for  the  deteimination  of  the  vertical  stress. 
These  charts  are  common  and  are  found  in  most  text  books 


on  Soil  Mechanics 
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A  graphical  solution,  which  proved  to  be  the  most  useful  was 
developed  by  Newmark  (19)*  This  method  provides  a  graphical  means  of 
computing  all  six  components  of  stress  under  any  loading  system  regular 
or  irregular  for  any  desired  value  of  Poisson1  s  ratio*  The  method  is 
based  on  the  Boussinesq  equations* 

Westergaard  Solution 

In  1938,  T/festergaard  (31)  published  a  solution  based  on  assumptions 
which  are  more  realistic  for  soil  than  those  of  Boussinesq*  He  assumes  a 
semi-infinite  material,  loaded  normal  to  the  surface*  The  material  is 
reinforced  with  numerous  closely  spaced  horizontal  sheets  which  retard 
lateral  motion*  His  solution  is  valid  for  any  Poisson* s  ratio,  and  is 
left  in  an  integral  form  for  expansion  to  suit  the  boundary  conditions 
of  a  particular  case. 

Taylor  (2U)  and  Fadum  (6)  have  simplified  the  Westergaard  analysis 
to  provide  a  quick  means  of  computing  the  vertical  stress  under  simple 
loading  conditions  for  Poisson* s  ratio  -  0.  This  is  done  by  means  of 
tables  and  charts  which  are  similar  to  those  published  by  these  authors 
for  solving  Boussinesq  equations.  The  vertical  stress  computed  by  these 
charts  are  approximately  two  thirds  of  the  value  computed  by  the 
Boussinesq  equations* 

There  seems  to  be  no  simplified  solutions  for  any  but  the 
vertical  stress  by  the  Westergaard  analysis. 
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Burmister  Solution 

The  solution  for  the  stress  due  to  a  load  applied  on  the  surface 
of  a  layered  system  was  first  published  by  Burmister  in  19k3  (2)*  He 
assumes  a  two  layered  system,  each  layer  being  homogeneous  and  isotropic 
but  with  different  physical  properties*  His  solution  shows  that  the  top 
layer  acts  as  reinforcing  and  distributes  the  load  over  a  wider  area 
onto  the  lower  layer*  Thus  the  stresses  in  the  lower  layer  computed  by 
this  method  are  lower  than  if  computed  by  Boussinesq*  Prior  to  this 
time  Casagrande  (3)  had  suggested  that  a  sand  and  gravel  layer  above  a 
softer  clay  layer  would  spread  the  load  and  distribute  the  stress  over 
a  wider  area  which  would  aid  in  smoothing  out  differential  settlements* 

The  Equations  Used 

The  consolidating  soil  at  the  Kitimat  Smelter  is  definitely  a 
layered  soil  and  for  this  reason  the  Westergaard  or  Burmister  solution 
would  probably  be  the  most  realistic*  The  compressible  material  also 
rests  on  solid  rock  at  a  depth  of  between  25  and  30  percent  of  the  width 
of  the  applied  load.  This  fact  would  suggest  an  approach  similar  to 
that  of  Pickett  and  others  who  have  considered  the  effect  of  a  rigid 
boundary  beneath  the  compressible  material*  None  of  these  methods  were 
chosen  however,  because  it  was  felt  that  the  accuracy  of  the  results  did 
not  warrant  the  time  and  effort  to  reduce  the  general  equations  into  a 
form  which  could  be  applied  to  this  particular  case*  It  was  therefore 
decided  to  use  the  Boussinesq  solution* 
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Of  the  many  available  approaches  using  the  Boussinesq  method  the 


following  were  considered  and  rejected# 

1*  Dividing  the  area  into  squares  small  enough  to  act  as 
point  loads  and  sum  the  stress  produced  by  each  "point 
loaded  square" •  It  can  be  easily  shown  that  to  be  within 
an  error  of  10$  the  distance  below  a  loaded  square  must  be 
approximately  three  times  the  width  of  the  square  before  the 
square  can  be  replaced  by  an  equivalent  point  load.  Stresses 
were  desired  within  the  clay  stratum,  from  about  15>0  to  300 
feet  below  the  applied  load.  Choosing  the  greater  depth, 
the  side  of  each  square  would  have  to  be  less  than  100 
feet,  and  the  number  of  squares  would  be  approximately  200. 

At  a  depth  of  15>0  feet  the  number  of  point  loaded  squares 
would  approach  1000.  This  number  could  be  reduced  somewhat 
because  the  contribution  of  the  outer  squares  to  the  total 
stress  would  be  insignificant  compared  to  that  of  squares 
close  to  the  point.  It  was  felt  however,  that  the  required 
number  of  computations  by  this  method  was  too  large  for 
practical  consideration. 

2.  The  use  of  JurgensGnrs  or  Fadum’s  tables.  These  tables  pro¬ 
vided  only  the  stress  components  in  one  plane  due  to  a  line 
or  symmetrical  load.  They  were  rejected  because  a  method 
to  apply  them  to  this  case  of  an  unsymme trie al  load  could 


not  be  found 
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3*  Integrate  the  Boussinesq  equations  as  published  by  Way  (30) 
over  the  boundary  of  the  applied  load  and  obtain  by  mathe¬ 
matical  computation  a  rigorous  solution  to  the  equations  at 
each  desired  point.  By  approximating  the  actual  applied 
load  by  a  wedge  the  equations  were  set  up  in  integral  form, 
but  their  complexity  discouraged  any  further  pursuit  along 
this  line. 

U*  The  use  of  Newmark’s  Influence  Charts.  As  previously  stated, 
these  charts  provide  a  simple  graphical  method  of  obtaining 
the  six  necessary  components  of  stress  due  to  any  normally 
applied  load,  and  for  any  value  of  Poisson* s  ratio.  As  this 
was  the  most  straight  forward  of  all  methods,  and  one  which 
was  most  readily  adaptable  to  this  particular  case  it  was 
adopted* 

Calculations 

Newmarks  influence  charts  are  graphical  solutions  to  the  Boussinesq 
equations.  There  is  a  separate  chart  for  each  of  the  six  components  of 
stress  (three  normal  and  three  shearing).  These  six  charts  are  set  up 
for  Poisson*s  ratio  of  0.5.  Three  other  charts  are  provided  to  correct 
the  applicable  stress  component  to  any  other  value  of  Poisson* s  ratio. 

The  use  of  the  charts  involves  drawing  the  loaded  area  on  tracing  paper 
to  an  appropriate  scale  and  superimposing  it  over  the  chart.  The  load 
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of  the  area  covering  each  square  is  multiplied  by  an  influence  value 
and  the  total  stress  is  the  summation  of  all  the  individual  squares. 

Having  determined  the  six  stress  components  from  Newmarks  charts 
the  direction  and  magnitude  of  the  major  principal  stress  can  then  be 
found  by  standard  equations  of  continuity  for  three  dimensional  stress 
system.  These  equations  are  derived  in  textbooks  on  the  theory  of 
elasticity  (for  example  Reference  28)  and  are  summarized  in  Figure  5U* 
After  the  direction  of  the  major  principal  stress  has  been  detemined 
its  horizontal  projection  in  the  direction  of  the  coordinate  axes  can 
be  easily  found  by  the  method  outlined  in  Figure  5^*  This  expresses 
tlie  horizontal  projection  of  major  principal  stress  as  a  fraction  of 
the  vertical  stress® 

The  vertical  settlement  which  has  taken  place  over  any  period  of 
time  was  accurately  known.  From  this  movement  the  horizontal  components 
of  the  maximum  movement  were  calculated.  This  involves  the  assumption 
tli at  the  vertical  and  maximum  movements  were  parallel  to  and  proportional 
to  the  vertical  and  major  principal  stresses.  The  consolidation  charac¬ 
teristics  in  the  horizontal  direction  were  probably  considerably  different 
from  those  in  the  vertical  direction,  but  lacking  information  on  this 
subject  the  soil  was  assumed  to  have  like  properties  in  all  directions. 

As  discussed  in  Chapter  V  a  detailed  set  of  observations  was 
made  for  column  lines  G  and  H  of  Potroom  2a.  Settlement  gauges  6,  7 
and  8  being  adjacent  to  column  line  H,  were  chosen  as  the  sites  for  the 
determination  of  theoretical  horizontal  stresses  and  displacements. 
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Stresses  at  depths  of  5>0,  15 0  and  250  feet  below  these  gauges  were 
obtained  for  Poisson* s  ratio  of  0.5*  Numerical  values  for  the  stress 
components  and  resulting  horizontal  displacements  are  shown  in 
Table  VIII*  To  show  the  variation  of  the  individual  stress  components 
with  depth  the  noimal  and  shearing  stresses  are  plotted  on  Figure  55» 
The  theoretical  horizontal  movement  expressed  as  a  percent  of  the 
vertical  movement  is  plotted  on  Figure  56.  From  the  total  building 
settlement  the  horizontal  movement  was  calculated  and  is  shown  on 
Figure  57* 

It  was  shown  in  Chapter  II  that  the  settlement  in  the  gravel 
layer  was  almost  completed  before  the  buildings  were  construe  ted# 
Settlement  at  any  point  within  the  gravel  layer  was  assumed  to  be  the 
same  as  the  surface  settlement.  In  Chapter  II  it  was  also  shown  that 
the  settlement  of  the  silt  and  clay  is  quite  similar.  The  proportion 
of  surface  settlement  assigned  to  any  point  within  the  silt  and  clay 
was  assumed  to  be  the  same  as  the  ratio  of  the  depth  below  the  top  of 
this  layer  to  the  total  thickness  of  these  compressible  soils. 

An  examination  of  the  individual  stress  components  shows  that 
the  major  principal  stress  is  only  slightly  greater  than  the  vertical 
stress.  The  shearing  stresses  are  very  small  and  the  direction  of  the 
major  principal  stress  is  close  to  being  vertical.  It  is  extremely 
difficult  to  make  exact  checks  on  the  computations  when  working  with 
three  figure  accuracy  and  small  magnitude  stresses  and  stress 
differences.  Mien  the  questionable  validity  of  assumptions  involving 
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the  use  of  the  Boussinesq  equations  for  this  soil  condition,  and  the 
method  of  graphical  computation  are  considered,  the  required  accuracy 
in  calculations  does  not  warrant  the  use  of  calculating  machines  and 
several  places  of  decimal*  She  resulting  curves  of  horizontal  displace¬ 
ment  are  plotted  as  they  were  calculated*  They  show  distinct  breaks 
which  are  probably  due  to  errors  incurred  from  calculations  involving 
differences  of  small  numbers*  This  in  no  way  invalidates  the  results* 
Sufficient  computations  were  made  to  indicate  the  general  magnitude 
and  trend*  None  of  the  results  showed  significant  disagreement  from 
this  average* 

To  show  the  general  trend  of  horizontal  movement  at  various  other 
places  over  the  Potline  1  and  2  area  stresses  and  displacements  were 
computed  at  a  depth  of  150  feet  below  each  settlement  gauge  in  this 
area*  The  results  are  shown  plotted  as  a  vector  diagram  on 
Figures  58  and  59* 

All  of  the  aforementioned  calculations  were  made  assuming 
Poisson1  s  ratio  equal  to  0*5*  This  assumes  no  volume  change  occurs 
during  elastic  deformation  which  is  the  usual  assumption  for  a 
saturated  soil.  Poisson* s  ratio  does  not  appear  in  the  Boussinesq 
equations  for  vertical  normal  stress  or  for  shearing  stress  on 
vertical  planes*  In  most  settlement  computations  the  vertical  stress 
is  the  only  component  required.  Since  the  Boussinesq  equations  are 

often  used  for  this  computation  the  choice  of  a  value  for  Poisson* s 
ratio  for  soil  is  often  not  important. 
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The  Westergaard  analysis  contains  Poisson1  s  ratio  in  the  expres¬ 
sion  for  vertical  stress.  The  value  is  usually  chosen  to  be  zero  which 
represents  no  lateral  expansion  due  to  compression  (6  and  2i|).  This 
gives  a  maximum  value  for  vertical  normal  stress.  Although  no  applic¬ 
able  method  of  using  the  Westergaard  equations  was  available  it  was 
decided  to  recompute  the  movement  at  a  few  points  from  the  Boussinesq 
equations  for  a  Poisson* s  ratio  equal  to  zero.  The  solution  for 
Poisson* s  ratio  of  0.5  and  zero  would  represent  limiting  values  for  the 
Boussinesq  solution.  It  was  felt  that  the  solution  for  Poisson* s  ratio 
of  zero  would  be  an  approximation  to  the  Westergaard  analysis.  Stresses 
and  horizontal  displacements  at  a  depth  of  150  feet  were  computed  for 
Poisson* s  ratio  of  zero  at  the  site  of  settlement  gauges  6,  7*  8,  9, 

10  and  11.  The  results  are  shown  on  Table  VIII  and  Figures  56  and 
57*  The  horizontal  movement  for  Poisson* s  ratio  equal  to  zero  was 
approximately  one  half  the  value  confuted  for  Poisson’s  ratio  equal 
to  0.5. 

Discussion  of  Results 

At  settlement  gauges  6,  7  and  8  the  calculated  horizontal  com¬ 
ponent  of  movement  was  within  the  range  of  10  to  25  percent  of  the 
vertical  movement  for  Poisson’s  ratio  of  0.5*  Using  Poisson’s  ratio 
equal  to  zero  reduced  this  horizontal  movement  by  approximately  one 
half.  The  calculated  movements  at  other  points  over  the  smelter  area 
showed  the  same  general  trend,  (Figure  58).  The  direction  of  movement 
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in  every  case  was  away  from  the  center  of  the  loaded  area  and  the  magni¬ 
tude  of  movement  increased  towards  the  edges  of  the  load. 

Jurgensen  (12)  has  published  a  diagram  showing  the  direction  of 
the  principal  stresses  at  several  points  below  and  at  some  distance 
beyond  the  edge  of  a  strip  footing.  The  direction  of  principal  stresses 
as  shown  herein  agreed  with  Jurgensen’s  diagram.  At  increasing  distance 
beyond  the  edges  of  the  load  Jurgensen fs  chart  shows  that  the  direction 
of  the  major  principal  stress  becomes  almost  horizontal. 

The  actual  horizontal  movement  at  each  point  was  computed  on 
the  basis  of  observed  total  settlement  since  construction  of  the 
buildings.  The  gravel  layer  was  considered  to  contribute  nothing  to  the 
vertical  settlement.  Thus  all  points  lying  within  the  depth  of  the 
gravel  layer  were  assumed  to  settle  the  same  as  at  the  surface.  Vertical 
settlement  of  points  within  the  compressible  silt  and  clay  layers  were 
reduced  in  proportion  to  their  depth  and  the  horizontal  movement  of 
these  points  was  computed  on  the  basis  of  this  reduced  vertical  settle¬ 
ment.  !&thin  the  depth  of  compressible  material  the  magnitude  of 
actual  movements  decrease  from  a  maximum  at  the  top  of  the  silt  to 
zero  at  the  rock  interface  because  zero  vertical  movement  was  assumed 
for  this  depth. 

The  magnitudes  of  horizontal  movement  in  the  upper  layer  ranged 
up  to  0*5  feet  at  the  edges  of  the  buildings  using  Poisson*s  ratio 
equal  to  0.5.  This  is  true  for  gauges  6,  7  and  8,  (Figure  57)  as  well 
as  the  other  gauges  over  the  smelter  area  (Figure  5 9)*  If  Poisson!s 
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ratio  equal  to  zero  were  used  it  would  be  expected  to  reduce  the  computed 
horizontal  movement  by  one  half,  or  to  approximately  0*25  feet  at  each 
end  of  the  buildings. 

Movements  at  the  upper  layer  were  based  on  calculation  of 
stresses  at  a  depth  of  £0  feet  below  the  bottom  of  the  fill#  This  is 
only  l/60th  of  the  width  of  the  loaded  area.  Theoretically  the 
Boussinesq  equations  hold  true  for  every  point  below  the  load  except  at 
the  contact  face.  It  can  easily  be  shown  that  for  simple  loading  con¬ 
ditions  the  induced  stress  changes  rapidly  toward  the  contact  face. 

While  the  equations  are  theoretically  correct  it  would  be  expected  that 
at  such  a  relatively  shallow  depth  great  care  and  precision  would  be 
required  to  obtain  accurate  results.  As  previously  mentioned  hoover, 
the  general  agreement  of  all  points  between  depths  of  50  and  250  feet 
would  indicate  that  the  results  are  sufficiently  reliable. 

The  actual  field  conditions  were  in  some  respects  considerably 
different  from  the  basic  assumptions  inherent  in  the  equations  which 
were  used  to  compute  the  stress  below  the  fill. 

1.  The  soil  was  not  an  isotropic  and  homogenous  material. 

There  were  three  principal  layers  of  quite  different 
elastic  properties,  and  within  each  layer  were  numerous 
horizontal  stratifications. 

.  The  rock  surface  at  a  relatively  shallow  depth  below  the 
load  presented  a  discontinuity  not  assumed  by  Boussinesq. 
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3o  The  main  movement  within  the  soil  was  due  to  consolidation 
rather  than  elastic  defoliation,  and  theories  of  elasticity 
may  not  strictly  apply. 

Modified  equations  have  been  introduced  by  Westergaard  and 
Bumister  to  attempt  to  account  for  horizontal,  stratification.  labile  no 
computations  were  made  for  this  particular  case  other  investigators  have 
suggested  that  these  horizontal  layers  tend  to  reduce  the  stress  within 
the  soil  below  the  value  computed  on  the  basis  of  a  homogeneous  and 
isotropic  material. 

Biot  and  Cummings  have  attempted  to  account  for  the  effect  of  a 
rigid  lower  boundary.  Cummings  suggested  that  the  stress  at  this 
boundary  could  be  approximated  by  using  Boussinesq  equations  and  a 
modified  depth  of  0.75  times  the  actual  depth.  He  also  suggested  that 
by  St.  Venantte  principle  the  stress  within  the  soil  a  short  distance 
from  the  rigid  boundary  would  be  the  same  as  computed  assuming  no  dis¬ 
continuity.  The  plots  of  stress  and  horizontal  movements  versus  depth 
show  little  variation  over  the  depth  considered.  Therefore  using  a 
somewhat  smaller  depth  would  have  little  effect  on  the  computed 
stresses. 

Conclusions 

The  direction  and  magnitude  of  the  major  principal  stress  was 
found  at  several  places  below  the  smeltersite.  The  direction  of  this 
stress  was  found  to  be  away  from  the  center  of  the  loaded  area.  It  was 
approximately  vertical  at  the  center  of  the  area  and  had  an  increasing 
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horizontal  inclination  toward  the  edges  of  the  fill.  The  direction  of 
the  major  principal  stress  did  not  vary  significantly  with  depth  within 
the  compressible  soil. 

By  assuming  that  movement  in  any  direction  due  to  consolidation 
was  proportional  to  the  component  of  the  major  principal  stress  in 
that  direction  horizontal  movements  were  calculated  at  various  places 
over  the  smeltersite.  The  vertical  component  of  movement  was  taken  as 
the  vertical  settlement  of  the  buildings  since  they  were  constructed. 
This  computed  horizontal  movement  was  away  from  the  center  of  the 
loaded  area.  At  the  surface  of  the  compressible  soils  it  varied  from 
zero  at  the  middle  of  the  area  to  a  maximum  of  0.5  feet  at  the  edges 
of  the  pot  room  buildings. 

It  was  shown  in  Chapter  II  that  settlement  since  construction 
of  the  buildings  has  been  almost  entirely  within  the  lower  silt  and 
clay  layers.  Therefore  this  computed  maximum  horizontal  movement 
would  occur  at  the  gravel  silt  interface,  which  is  50  to  120  feet  below 
the  surface.  How  this  outward  movement  is  transmitted  through  the 
gravel  layer  is  uncertain.  Certain  physical  boundary  conditions 
must  be  met.  Some  of  these  ares 

1.  A  closed  system  -  no  cracks  can  form  within  the  soil 
at  depth. 

2.  The  horizontal  movement  must  be  zero  at  the  rock  surface. 
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3»  Since  there  was  no  significant  settlement  within  the  upper 
gravel  layer,  approximately  5>0  to  120  feet  thick,  the  seat 
of  all  horizontal  movement  was  within  the  silt  and  clay 
below* 

Ij.*  There  cannot  be  a  discontinuity  at  the  silt-gravel  interface 
and  the  movement  at  each  layer  must  be  the  same  in  magnitude 
and  direction* 

In  Chapter  V  it  will  be  shown  that  the  ends  of  the  potroom 
buildings  appear  to  have  moved  closer  together  by  approximately  one  half 
foot*  This  would  represent  a  movement  of  J  foot  at  each  end  of  the 
potroom,  or  approximately  one  half  of  the  movement  computed  on  the 
basis  of  a  non  vertical  major  principal  stress*  The  direction  of  the 
observed  movement  however  was  inward,  and  opposite  to  the  outward 
movement  calculated  from  the  stress  consideration.  An  explanation  for 
this  apparent  contradiction  is  presented  in  Chapter  V  along  with  a 
discussion  of  other  factors  related  to  the  problem  of  horizontal 
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CHAPTER  IV 


MODEL  STUDY 

Introduction 

To  tiy  to  obtain  a  general  picture  of  the  behaviour  of  the  soil 
surface  as  a  result  of  settlement  it  was  decided  to  construct  a  model# 
This  was  done  by  substituting  elastic  bands  for  the  compressible  soil, 
and  inducing  in  them  a  tensil  rather  than  a  compressive  stress#  This 
would  give  results  exactly  opposite  to  the  movement  in  the  field# 

Description  of  Model 

The  model  study  was  made  prior  to  receiving  detailed  information 
regarding  the  movements  along  column  lines  G  and  H.  Column  line  A 
was  chosen  for  the  study#  Along  this  column  line  the  subsoil  and 
settlement  profiles  were  plotted#  The  settlement  profile  was  plotted 
above  the  subsoil  with  settlements  increasing  •upwards  as  shown  in 
Figure  60#  The  same  horizontal  scale  was  used  in  each  case,  but  to 
accentuate  the  results  vertical  scales  were  distorted#  These  scales 
were: 

Subsoil:  Horizontal  lw  s  5 0  feet. 

Vertical  ltJ  -  30  feet# 

Settlement  Profile:  Horizontal  ln  ■  5>0  feet. 


Vertical  1"  s  1  foot 
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Since  only  qualitative  results  were  expected  from  the  test,  great 
precision  was  not  required  in  setting  up  the  apparatus#  For  this  reason 
the  rubber  bands  which  substituted  for  the  soil  were  chosen  on  the  basis 
of  availability#  Three  sizes  were  used  for  each  of  the  three  subsoil 
layers# 

\  inch  wide  for  gravel 
3/l6  inch  wide  for  silt 
l/8  inch  wide  for  clay. 

.All  the  bands  were  approximately  l/ 2l|th  inch  thick#  Within  the  range 
of  strain  applied  the  ratio  of  the  strains  was  constant  and  equal  to 
J"  (Gravel):  3/16”  (Silt);  l/8»  (Clay)  =1  :  1.2  :  1.? 

Using  the  appropriate  size  band  in  each  soil  layer  a  grid  of  1.8 
inch  squares  was  formed  with  horizontal  and  vertical  orientation.  The 
vertical  side  of  each  square  fell  immediately  below  the  column  on  which 
settlement  observations  were  taken.*  The  arrangement  was  such  that  the 
soil  surface  was  replaced  by  a  horizontal  \  inch  band  but  the  lower 
squares  were  sometimes  not  complete#  Seven  inch  long  wire  connectors 
stretched  from  the  line  of  zero  settlement  on  the  settlement  pattern, 
to  the  line  of  zero  deformation  on  the  soil  profile  where  they  were 
joined  to  the  rubber  bands.  The  horizontal  bands  extended  a  variable 
distance  beyond  the  end  of  the  building.  Figure  60  shows  these  dis¬ 
tances  to  the  same  horizontal  scale  as  the  model.  The  entire  model 


*  See  figure  3  for  location  of  settlement  observation  points. 
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was  constructed  on  a  sheet  of  2^*  x  3*  x  \  inch  plywood*  The  rubber 
bands  were  given  an  initial  strain  equal  to  20  percent  of  their  original 
length*  They  were  then  pinned  to  the  boundaries  as  described.  The 
corners  of  the  squares  were  then  attached  by  rubber  canent* 

Operation  of  the  Model 

When  the  cement  had  hardened  the  wire  connectors  were  moved  from 
the  line  of  zero  settlement  to  the  line  showing  the  total  settlement  of 
the  building,  and  pinned  in  place*  Thus  the  soil  surface  at  each 
settlement  observation  point  was  given  an  induced  vertical  displace¬ 
ment  corresponding  to  the  actual  conditions  without  introducing  any 
horizontal  fixity* 

When  the  soil  had  been  deformed  as  described  the  horizontal 
movement  at  the  ground  surface  was  observed  at  every  settlement 
observation  point.  This  process  was  repeated  with  the  horizontal  bands 
pinned  at  various  distances  beyond  the  edges  of  the  building* 

Observations 

The  results  of  the  model  study  tests  are  summarized  on  Figure  60* 
The  horizontal  movements  shorn  were  measured  to  the  scale  of  the  model, 
but  were  really  only  a  matter  of  one  tenth  inch  and  less.  It  can  be 
seen  that  the  movements  which  occurred  in  the  model  indicated  the  end 
column  footings  of  the  prototype  moved  together.  The  greatest  movement 
was  in  the  vicinity  of  maximum  settlement. 
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Conclusion 

The  model  was  not  constructed  with  great  precision,  and  the 
results  can  only  be  considered  of  a  qualitative  nature.  The  results 
however  do  indicate  a  horizontal  movement  of  the  ground  surface*  This 
movement  was  generally  toward  the  center  of  the  potroom  buildings. 
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CHAPTER  V 


GEOMETRY  OF  SUBSIDENCE 

Description  of  Measurements 

The  first  indication  of  serious  structural  damage  in  the  Potroom 
buildings  -which  could  be  attributed  to  settlement  came  in  the  summer  of 
1957*  Many  of  the  steel  columns  had  developed  a  distinct  tilt.  It  was 
first  thought  that  the  tendency  was  for  the  tops  of  the  columns  to  tilt 
towards  the  west,  and  it  seemed  probable  that  they  were  being  given  a 
horizontal  thrust  in  this  direction  at  the  height  of  the  crane  rail. 

The  most  serious  effect  was  on  the  west  one  half  of  the  buildings,  but 
there  was  some  evidence  of  an  eastward  thrust  of  the  column  tops  on  the 
east  end  as  well.  Associated  with  this  apparent  movement  of  the  steel 
superstructure  were  buckled  diagonal  bracing  members,  and  cracked  con¬ 
crete  supports.  Longitudinal  and  cross-sections  of  a  typical  Potroom 
are  shown  in  Figures  5  and  6. 

•Since  the  cause  of  the  damage  was  not  manifest  by  general  obser¬ 
vations  and  measurements  made  at  the  time,  a  detailed  set  of  measure¬ 
ments  were  taken  during  November  and  December  1957*  Column  lines  G  and 
H  (Potroom  2A)  were  chosen  for  the  detailed  study  because  they  were 
situated  in  an  area  of  maximum  differential  settlement.  The  following 
measurements  were  taken  on  each  column  line. 
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1*  Horizontal  distance  between  Columns  1  and  lj.8  at  the  elevation 
of  the  steel  base  plate*  (i*e.  top  of  the  concrete  columns) 
and  at  the  elevation  of  the  working  floor* 

2.  Horizontal  distance  between  the  center  line  of  each  column 
at  the  elevation  of  the  steel  base  plate* 

3*  Deviation  from  plumb  taken: 

(i)  Between  the  top  of  the  concrete  column  and  7.6  feet 
beloxtf* 

(ii)  Between  steel  base  plate  and  1  inch  below  the 
gusset  plate  on  the  crane  beam  web. 

(iii)  Between  the  steel  base  plate  and  the  top  of  the 
steel  column. 

ii*  Elevation  of  the  structural  steel  base  plates. 

5>*  Distance  across  expansion  joints  at  the  crane  beam. 

The  end  to  end  distances  were  each  measured  twice  at  the  eleva¬ 
tion  of  the  base  plates  and  once  at  the  working  floor*  A  further 
check  was  made  by  aligning  the  end  columns  by  transit  and  noting 
the  difference  in  length*  A  steel  chain,  calibrated  to  a  standard 
chain  was  used  and  corrections  for  sag  tension  and  temperature  were 
applied.  All  measurements  were  taken  inside  the  potroom  where  the 
radiant  heat  from  the  smoldering  pots  created  variable  temperatures 
along  the  room. 

The  first  measurement  was  made  November  25,  1957  at  the  eleva¬ 
tion  of  the  steel  base  plates*  The  temperature  was  recorded  at  the 
23  and  73  foot  marks  of  the  100  foot  steel  chain*  On  November  26 
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measurement  was  taken  at  the  elevation  of  the  working  floor  with 
temperatures  similarly  recorded*  The  end  points  of  this  measurement 
were  vertically  below  that  of  the  previous  day,  and  the  distance  was 
between  the  center  lines  of  the  end  steel  columns* 

On  December  5,  a  further  measurement  was  made  at  the  elevation 
of  the  steel  base  plates  and  air  temperatures  were  recorded  at  every 
column*  These  temperatures  varied  from  1*0  to  81*°  F  with  as  much  as 
10°  F  difference  between  adjacent  columns*  A  comparison  of  all 
measurements  is  shown  in  Table  IX*  The  difference  in  length  as 
measured  by  the  transit  alignment  was  assumed  to  be  correct.  On  this 
basis  the  probable  error  in  measuring  the  overall  length  varied  from 
.060  feet  on  November  2$  to  *00?  feet  on  November  26* 

Measurement  of  the  lean  of  the  concrete  columns  was  done  by 
suspending  a  lead  plumb  bob  from  a  wooden  frame  held  to  the  side  of 
the  column.  Allowance  was  made  for  any  deviation  in  width  of  the 
column.  Measurements  were  taken  to  the  nearest  l/l6  inch  and  are 
probably  accurate  to  the  nearest  l/8  inch.  Lead  and  wood  were  chosen 
because  they  would  not  be  effected  by  the  extremely  high  electro -magnetic 
forces  in  the  potroom. 

Measurement  of  the  tilt  of  the  steel  column  was  done  by  means 
of  a  transit  and  reading  the  intersection  of  the  vertical  cross  hair 
on  a  steel  tape  held  against  the  base  of  the  column*  Again  measurements 
are  probably  accurate  to  the  nearest  l/8  inch. 
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The  elevation  of  the  steel  base  plates  were  compared  to  the 
recommended  setting  elevation.  In  this  way  the  total  settlement  of 
each  column  since  construction  was  determined.  The  elevations  are 
probably  accurate  to  the  nearest  0.01  feet  but  since  no  previous 
check  was  made  of  the  exact  setting  elevation,  and  as  all  columns  may 
not  have  been  brought  to  their  final  position  the  day  the  recommended 
setting  was  given  it  is  probable  that  some  of  the  settlement  values 
may  be  in  error  by  a  maximum  of  0.05  feet. 

The  measurement  across  the  expansion  joints  were  mainly  for  the 
purpose  of  future  comparisons  and  gave  only  an  indication  of  whether 
the  expansion  was  fully  taken  up  or  not.  Inttiile  an  indication  of  the 
relative  amount  the  joint  is  open  or  closed  could  be  observed,  it  was 
not  possible  to  tell  by  visual  inspection  whether  the  joints  were  fully 
closed.  A  power  failure  which  occurred  December  1957  resulted  in  a 
lowering  of  the  air  temperatures  inside  the  pot  rooms  by  an  estimated 
10  to  15°  F.  Measurement  of  the  expansion  joints  before  and  after 
showed  no  significant  change  except  at  G30  and  H30  which  opened  3/8 
and  7/16  inch  respectively.  A  few  measurements  of  the  lean  of  steel 
columns  before  and  after  the  power  failure  showed  movement  of  up  to 
J  inch  but  insufficient  readings  were  taken  to  establish  a  pattern. 

The  results  of  all  measurements  are  summarized  on  two  large  drawings. 
Figures  68  and  69  enclosed  herewith. 
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Interpretation  of  Measurements 

The  most  significant  fact  emerging  from  the  detailed  measurements 
is  that  both  column  lines  G  and  H  appeared  to  have  shortened  approxi¬ 
mately  one  half  foot*  Two  explanations  seemed  possible* 

1*  The  columns  were  initially  set  in  error  by  this  amount. 

2*  The  settlement  has  drawn  the  columns  closer  together. 

Only  a  few  checks,  independent  of  the  contractors  own  work  were 
made  at  the  time  of  construction.  These  checks  however  indicated  that 
the  initial  column  setting  was  correct.  An  error  of  one  half  foot  in 
1100  feet  (approximately  .Oi;  feet  in  a  100  foot  chain)  is  rather  large* 
It  is  improbable  that  an  error  in  the  chain  of  this  magnitude  would  go 
unnoticed  at  the  time  of  construction* 

Other  factors  also  indicate  that  the  initial  setting  was  correct, 
and  the  column  bases  have  since  moved  to  their  present  position*  lniien 
the  power  was  shut  off  and  the  temperatures  lowered  app roximately  15° 
the  1100  foot  crane  beam  would  have  tended  to  shorten  app  roximately 
0.10  feet**'  The  fact  that  only  the  joint  at  column  30  showed  any 
significant  opening  would  indicate  that  a  compressive  force  existed 
across  the  other  joints  which  was  not  relieved  by  the  temperature 
change.  Movement  of  the  steel  column  tops  during  this  temperature 

*  Coefficient  of  linear  expansion  of  steel  s  6.7  x  10~°  ft./ft./°F, 

Steel  Construction,  AISC,  New  York,  19\±9,  p*3U8. 
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change  may*  also  indicate  that  a  stress  existed  in  the  crane  beam  and 
was  transferred  across  the  expansion  joints.  An  examination  of  the  gen¬ 
eral  lean  of  the  steel  columns  showed  that  while  the  leans  were  very 
erratic*  it  appeared  that  the  tops  were  being  pushed  away  from  the 
zone  of  maximum  settlement.  This  movement  of  the  steel  columns  and 
compression  of  the  crane  beam  would  probably  only  occur  if  the  base  of 
the  columns  moved  together. 

The  horizontal  thrust  transferred  to  the  crane  beam  when  the 
crane  stops  and  starts  at  the  ends  of  the  potroom  has  also  been  sug¬ 
gested  as  a  possible  explanation  for  some  of  the  lateral  movement  of 
the  steel  column  tops.  "While  this  may  be  a  contributing  factor  to 
the  distortion  of  the  superstructure  there  are  other  more  probable 
causes.  A  factor  which  will  be  discussed  later  in  greater  detail  is 
that  the  major  shortening  and  the  greatest  amount  of  column  disturbance 
lies  in  the  region  of  maximum  differential  settlement.  This  suggests  a 
relationship  between  settlement,  column  distortion  and  horizontal 
movement. 

From  the  above  reasoning  it  appeared  likely  that  the  deviation 
from  design  spacing,  and  the  distortion  of  the  superstructure  were 
related  in  part  at  least  to  settlement  of  the  buildings.  The  problem 
was  therefore  to  find  an  explanation  for,  or  relationship  between 
settlement,  horizontal  movement,  and  column  disturbance.  A  theory  is 
suggested  to  account  for  the  horizontal  movement  of  the  ground  surface, 
but  a  detailed  study  of  the  distortion  of  the  superstructure  was  not 
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Surface  Geometry 

The  first  attempt  to  explain  the  shortening  of  the  column  line  was 
to  consider  the  change  in  length  due  to  simple  geometry*  The  plot  of 
total  settlement  of  each  column  hereafter  referred  to  as  the  subsid¬ 
ence  curve,  was  approximated  by  two  straight  lines  joining  each  end 
to  the  point  of  maximum  settlement*  The  two  triangles  fomed  by  these 
sloping  lines,  the  common  vertical  and  the  original  horizontal  surface 
were  considered.  By  the  Pythagoras  theorem  the  length  of  the  horizontal 
side  was  computed  assuming  the  sloping  length  remained  constant  during 
the  settlement*  It  was  found  by  this  method  that  the  combined  shorten¬ 
ing  would  be  less  than  *005  feet.  Obviously  another  explanation  was 
required. 

Simil ar  Problems  from  other  Sources 

Literature  on  the  subject  of  settlement  and  subsidence  of  the 
ground  contains  many  examples  of  similar  movements  and  resulting 
structural  damage. 

For  many  years  the  mining  industry  has  been  concerned  because  of 
subsidence  at  the  surface  over  a  mining  area.  The  concern  in  this  case 
is  generally  because  many  mines  go  beneath  important  existing  structures 
which  are  susceptible  to  damage  when  the  earth  subsides.  It  has  been 
observed  that  horizontal  movements  of  the  earth  !s  surface  is  usually 
associated  with  vertical  subsidence,  and  it  is  these  horizontal  move¬ 
ments  which  are  most  damaging  to  existing  structures.  Unfortunately 
much  of  the  literature  on  this  subject  was  in  European  publications. 
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mostly  English  and  German,  and  not  generally  available  in  this  country. 
The  reason  was  perhaps  that  with  the  European  population  being  more  dense 
than  in  America,  especially  over  the  mining  area,  the  problem  is  more 
acute  in  that  part  of  the  world.  Much  of  the  literature  contained  only 
a  general  description  of  ground  subsidence  at  a  particular  area.  How¬ 
ever,  many  theories  and  explanations  have  been  advanced  to  attempt  to 
explain  the  horizontal  movements  which  often  occur  when  the  ground 
subsides  (20).  The  fact  that  new  theories  are  still  being  advanced 
would  indicate  that  much  is  yet  to  be  learned.  A  thorough  understand¬ 
ing  seemed  to  be  hampered  somewhat  because  of  lack  of  reliable  data. 
Almost  all  writers  stressed  the  need  of  more  precise  measurements, 
both  horizontal  and  vertical,  over  a  subsiding  area. 

Vertical  and  horizontal  movements  of  the  ground  surface  have  also 
been  a  serious  problem  in  areas  where  the  water  table  has  been  lowered, 
or  where  oil  has  been  pumped  from  its  reservoir  in  the  sand  beneath  the 
surface.  In  fact  one  author  states  that  he  knows  of  no‘ older  oil  field 
where  sufficient  data  will  not  show  subsidence^1  (23) • 

To  indicate  the  similarity  of  the  problem  at  Kitimat  with  other 
subsiding  areas  a  few  examples  will  be  cited.  In  a  recent  textbook 
Henry  (10)  has  noted  several  cases  of  ground  subsidence  with  resulting 
horizontal  movement.  In  the  Netherlands  a  building  1+20*  long  with 
columns  every  1*0 1 -10”  subsided  approximately  II  feet  and  had  a  maximum 

differential  settlement  of  6  inches  between  the  ends  of  the  building. 

The  maximum  horizontal  movement  was  3  inches  per  column  bay.  \$ien  the 
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differential  settlement  was  lj-  inches  over  the  length  of  the  building, 
horizontal  movements  of  plus  and  minus  one  inch  per  bay  were  observed* 
The  seat  of  settlement  was  the  mining  of  a  nearby  coal  seam  2f  to  6  feet 
thick  and  1300  to  2200  feet  below  the  surface*  Another  example  also 
cited  by  Henry  was  a  horizontal  movement  of  over  two  feet  due  to 
mining  a  coal  seam  approximately  2000  feet  below  the  surface. 

Gibson  (7)  has  reviewed  some  general  data  regarding  subsidence 
over  a  coal  mining  area.  He  gives  the  following  general  observations; 
depth  of  working  -  2£0  yards,  seam  thickness  -  3* -5 11 5  total  subsidence  - 
2!-9n,  maximum  slope  of  the  ground  surface  1.61  inches  in  10  feet  and 
horizontal  movements  of  the  order  of  1|  inches  in  100  feet. 

An  example  of  subsidence  due  to  lowering  the  ground  water  table 
has  occurred  at  Houston  Texas  (1?)  where  a  30  mile  diameter  area  has 
subsided  a  maximum  of  3»5>  feet.  This  has  been  accompanied  by  faulting 
and  horizontal  movements  of  from  6  to  16  inches. 

Shoemaker  and  Thor  ley  (23)  have  described  a  maximum  subsidence  of 
20  feet  over  the  Wilmington  oil  field  of  Long  Beach,  California*  The 
affected  area  is  approximately  5  x  10  miles  and  horizontal  movements 
have  been  of  the  order  of  2  to  b  feet.  Speaking  of  the  same  general 
area  Henry  cited  an  example  of  11  feet  vertical  subsidence  accompanied 
by  a  5’  shortening  of  a  6000  foot  base  line. 

The  most  serious  effects  of  ground  subsidence  appeared  to  be  a 

result  of  the  horizontal  movements  of  the  earths  surface  (7>  10  and  23). 
Structures  which  have  been  adversly  affected  are  buildings,  sewers  and 
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drains,  pavements,  bridges,  rivers,  canals  and  railroads. 

Henry  stated  that  for  buildings  at  least,  horizontal  movements 
can  be  discerned  before  vertical  movements  yet  they  and  their  associa¬ 
ted  horizontal  forces  are  more  difficult  to  account  for  in  the  structural 
design  than  vertical  movement.  References  7*  10  and  11  describe  various 
ways  of  designing  structures  to  withstand  subsidence.  In  general,  all 
were  agreed  that  the  most  stable  structure  for  a  subsiding  area  is  one 
which  is  free  to  move  with  the  ground.  A  large  number  of  closely 
spaced  expansion  joints  and  flexible  connections  are  recommended  to 
allow  for  considerable  horizontal  and  vertical  movement.  They  also 
recommend  that  between  the  joints  the  structure  be  as  rigid  as 
possible. 

Description  of  the  Subsidence  Curve 

The  available  literature  was  agreed  on  the  physical  characteristics 
of  ground  subsidence  caused  by  mining,  and  the  characteristics  of  sub¬ 
sidence  curves  from  other  causes  seemed  to  follow  this  same  pattern.  A 
good  description  of  the  subsidence  curve  is  given  by  Rellensmann  (22) 
and  is  summarized  on  Figure  61.  He  states  that  the  horizontal  shifting 
is  always  toward  the  region  of  greatest  subsidence,  and  the  maximum 
shift  occurs  at  the  point  of  greatest  slope  of  the  subsidence  curve. 

There  is  no  horizontal  movement  at  the  middle  of  the  subsidence  trough 
and  outside  the  zone  of  vertical  settlement.  This  results  in  tension  or 
stretching  of  the  soil  near  the  outside  edges  of  the  curve  and  compression 
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or  shortening  toward  the  zone  of  maximum  settlement*  Since  soil  can  take 
little  tension  the  outside  zone  often  contains  cracks*  Rellensmann 
points  out  however  th^t  the  cracks  may  be  too  small  to  be  seen  and 
their  presence  is  indicated  only  by  measurements* 

Theories  have  been  developed  to  attempt  to  calculate  in  advance 
the  magnitude  of  vertical  and  horizontal  movement*  Rellensmann  states 
that  most  of  this  work  has  been  done  in  Germany*  King  and  Smith  (13 ) 
performed  model  tests  to  estimate  surface  movements.  This  was  done  by 
removing  small  blocks  from  the  bottom  of  a  container  of  sand.  Their 
results  were  only  of  a  qualitative  nature  but  agreed  with  actual  observed 
cases.  In  a  recent  paper  Perz  (20)  has  summarized  the  attempts  to  solve 
the  subsidence  problems  under  the  following  headings. 

1.  Statistical  evaluation  of  measurements  —  good  for  only 
one  district  from  which  measurements  came. 

2*  Mathematical  treatment  based  on  certain  assumptions. 

3.  Precalculation  based  on  properties  of  overlying  rock  — 

limited  because  of  inadequate  knowledge  of  rock  properties. 
k*  Model  tests  —  questionable  if  they  can  reproduce  field 
conditions . 

He  suggests  that  although  there  has  been  numerous  attempts  to  solve 
the  problem  none  are  adequate.  He  proposes  a  solution  based  on 
"general  logic"  along  with  the  basic  assumptions  for  horizontal 

movement  as  shown  in  Figure  62.  He  assumes  the  soil  within  the  area 
of  influence  to  be  composed  of  a  number  of  triangular  wedges* 
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When  ore  is  removed  these  wedges  slide  into  the  void*  Because  the  area 
of  influence  is  wider  at  the  surface  than  at  the  ore  seam  the  sliding 
wedges  have  a  lateral  component  of  motion  as  shown*  By  assuming  that 
the  soil  can  withstand  some  tension  and  compression  thus  giving  zero 
motion  at  the  edges  and  at  the  middle  of  the  subsidence  trough  a  horiz¬ 
ontal  displacement  curve  similar  in  shape  to  that  of  Figure  61  is  obtain¬ 
ed*  The  expression  for  vertical  movement  contains  several  physical 
constants  to  be  evaluated  from  particular  cases  and  is  not  of  interest 
to  the  Kitimat  problem* 

The  equation  for  horizontal  movement  is  based  on  the  geometric 
shape  of  the  subsidence  curve  and  the  depth  to  the  mine  seam*  This 
procedure  of  fitting  certain  calculated  curves  to  observed  data 
appeared  to  be  the  method  of  solution  used  by  most  investigators* 

l/dhile  the  Perz  solution  was  interesting  it  was  not  expected 
that  the  results  from  its  use  would  give  perfect  agreement  to  the 
observed  data  at  Kitimat*  The  main  reason  was  because  the  solution  is 
based  on  removal  of  material  from  a  well  defined  rectangular  volume  and 
contains  as  a  major  factor  the  depth  h  from  the  surface  to  this  volume. 

The  seat  of  settlement  at  the  Kitimat  Smeltersite  was  at  some 
considerable  distance  below  the  surface  but  the  boundary  was  not  well 
defined,  and  the  volume  decrease  within  the  soil  was  not  a  simple 
rectangular  shape  as  it  is  with  the  Perz  solution  of  the  mining 
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Horizontal  Displacement  Caused  by  Soil  Consolidation 

Although  the  equations  based  on  mining  subsidence  may  not  strictly 
apply  to  the  problem  at  Kitimat,  there  are  certain  geometric  similarities 
which  suggest  that  the  problems  are  related.  The  seat  of  settlement  was 
at  an  appreciable  depth  below  the  surface,  and  the  subsidence  curves 
were  similar.  It  has  been  shown  that  very  little  settlement  has 
occurred  within  the  gravel  layer  since  the  buildings  were  constructed. 
This  layer  including  the  fill  is  from  5>0  to  120  feet  thick. 

Basic  As sumptions 

On  the  basis  of  the  slope  of  the  subsidence  curve  at  Kitimat  a 
possible  method  of  explaining  the  associated  horizontal  movements  is 
suggested.  The  basic  assumptions  are  as  follows  and  the  physical 
conditions  are  shown  on  Figure  63 • 

1.  No  settlement  occurs  in  the  upper  gravel  layer. 

2.  The  line  between  the  firm  gravel,  and  the  consolidating 
soil  beneath  is  well  defined,  and  at  depth  h  from  the 
surface. 

3*  The  vertical  movement  at  the  surface  is  equal  to  that  at 
depth  h,  hence  the  subsidence  curve  at  depth  h  is  parallel 
to  that  at  the  surface. 

1±.  The  gravel  can  withstand  sufficient  tension  and  compression 
strains  to  give  zero  lateral  motion  at  the  point  of  zero. 


and  maximum  settlement 
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$•  Perfect  friction  exists  at  the  gravel-silt  interface,  and 
no  sliding  can  occur  along  this  surface. 

60  No  cracks  can  form  within  the  soil  which  do  not  originate 
at  the  surface. 

These  assumptions  were  of  course  not  completely  met  in  the  field. 

A  small  amount  of  settlement  will  occur  within  the  upper  gravel  layer. 

The  boundary  between  the  gravel,  and  the  consolidating  soil  will  not  be 
well  defined.  As  explained  by  Taylor  (2k)  the  soil  near  the  drainage 
face  consolidates  at  a  faster  rate  than  the  soil  in  the  middle  of  the 
layer.  Thus  as  consolidation  proceeds  the  depth  h  will  increase. 
Assumptions  3  to  6  were  probably  in  close  agreement  with  field  conditions. 

Mechanism  of  Movements 

Consider  any  element  of  soil  ABCD  within  the  gravel  layer 

(Figure  63 )•  As  a  soft  soil  beneath  this  element  consolidates,  C  moves 

- 1 

vertically  down  to  C*  and  D  moves  to  Dh  Since  point  D  settles  more 
than  C  the  new  position  of  the  base  of  the  element,  CD’,  will  be 
inclined  at  an  angle  o(  to  the  original  base  CD.  Since  there  can  be  no 
slipping  or  tearing  at  the  base  of  this  element  ‘the  line  G’Df  must  be 
vertically  below  CD.  To  keep  the  angle  between  the  base  and  sides  of  the 
element  constant  the  surface  points  A  and  B  must  move  laterally  as  well 
as  vertically.  By  simple  geometry  it  follows  that  for  this  assumption 
to  be  fulfilled  the  horizontal  movement  v  of  any  point  can  be 
expressed  by 

(  7) 


v  =  h, 


? 


' 


: 


. 


1  1 


O  •  . 


85 


where  oC  is  the  angle  between  the  initial  and  final  position  of  the 
base  of  the  element*  Assuming  the  subsidence  at  the  surface  to  be 
parallel  to  that  at  the  base,  then  oC  is  the  angle  between  the  initial 
surface  and  the  subsidence  curve* 

Application  to  the  Kitimat  Problem 

The  following  assumptions  were  made  when  applying  the  above 

equation  to  the  conditions  at  Kitimat 

h  =  total  depth  of  fill  and  natural  gravel  as  shown  on 
Figure  2* 

ot  s  the  slope  of  the  curve  of  settlement  of  buildings 
since  their  construction. 

Computations  for  the  horizontal  movement  by  this  method  are  shown  on 
Tables  X  and  XI*  The  results  are  shown  plotted  on  Figures  6U  and  65* 
These  plots  also  show  a  comparison  to  the  actual  measurements* 

It  was  previously  stated  that  the  measurements  between  each 
column  was  made  only  at  the  elevation  of  the  steel  base  plate*  Figures 
68  and  69  were  used  to  determine  the  distance  between  the  columns  at  the 
ground  surface.  Each  concrete  column  was  plotted  in  its  relative  posi¬ 
tion  and  at  the  measured  lean.  The  departure  from  design  distance 
between  the  columns  at  the  ground  surface  was  then  measured  directly 
from  the  plot.  Some  starting  point,  which  could  be  assumed  to  have  not 
moved  laterally  was  required  in  order  to  compute  the  total  horizontal 
movement  of  any  particular  column.  According  to  all  the  available 
information  on  ground  subsidence  the  three  possible  choices  were  at  the 
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center  of  the  subsidence  trough  or  point  of  maximum  settlement,  and 
beyond  either  edge  of  the  subsidence  effect*  Having  no  measurements 
beyond  the  zone  of  settlement,  the  approximate  center  and  lowest  spot 
of  the  subsidence  trough  was  chosen  as  the  point  of  zero  lateral 
motion*  In  each  direction  from  this  point  the  deviations  from  design 
spacing  were  accumulated  to  give  the  total  lateral  movement  of  each 
column  base*  On  the  same  plot,  Figures  6k  and  65,  is  shown  the  lateral 
movement  as  computed  from  the  Perz  equation  for  mining  subsidence* 

These  results  show  several  interesting  and  significant  factors* 

1*  The  point  of  zero  lateral  motion  is  common  to  both  solu¬ 
tions  and  agrees  with  the  actual  measurements. 

2*  There  is  close  agreement  at  the  end  columns. 

3*  The  computed  movement  at  the  intermediate  columns  appears 
to  fluctuate  considerably* 
l±*  The  Perz  equation  gives  results  far  too  large. 

The  first  observation  was  expected  since  it  is  a  basic  assumption  of 
the  Perz  equation,  and  of  the  calculations  from  observed  data.  Since 
the  subsidence  curve  was  horizontal  at  this  approximate  location  the 
horizontal  movement  was  therefore  zero. 

The  most  significant  factor  is  the  close  agreement  of  the  end 
columns.  Since  the  overall  distance  was  more  accurately  known  than  the 
distances  between  individual  columns  the  close  agreement  indicated  the 

correct  nature  of  the  equation*  Fluctuation  of  movement  at  intermediate 
columns  is  probably  partly  due  to  inaccuracies  in  determining  the  actual 


subsidence  curve.  To  smooth  out  these  irregularities  the  results  of  a 
computation  based  on  average  major  changes  in  slope  is  shown  on 
Figures  66  and  6?*  These  curves  show  fair  agreement,  especially  in 
geometric  shape,  with  observed  values  at  intermediate  as  well  as  at 
the  end  columns.  It  should  also  be  noted  that  both  observed  and  computed 
curves  show  that  most  of  the  horizontal  movement  has  taken  place  in  the 
region  of  maximum  slope  of  the  subsidence  curve. 

Discussion  of  Results 

The  deviation  between  the  observed  and  calculated  values  may  be 
due  to  several  causes. 

1.  Inaccuracies  in  the  initial  column  setting,  both  vertically 
and  horizontally.  Small  deviations  from  the  recommended 
setting  position  would  go  unnoticed  in  the  field,  but  would 
induce  an  error  into  the  calculated  movement. 

2.  Inaccuracy  in  determining  the  actual  movement  at  the  ground 
surface.  These  values  were  computed  on  the  basis  of  hori¬ 
zontal  distances  measured  at  the  top  of  the  concrete 
column  and  the  slope  of  the  column  to  the  ground. 

3©  Incorrect  choice  of  the  depth  of  gravel  h.  This  value  was 
based  on  Figure  2  which  was  in  turn  based  on  a  comparatively 
few  widely  spaced  test  holes.  The  value  for  h  could  probably 
be  considered  as  somewhat  greater  than  the  depth  of  gravel 
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since  the  underlying  silt  layer  would  probably  be  nearly 

completely  consolidated  for  some  distance  from  this  surface. 

The  equation  does  account  for  any  tension  or  compression  forces  in 
the  gravel  due  to  rotation  of  these  elements  about  a  neutral  axis  at 
the  base.  It  would  seem  that  these  forces  would  in  effect  tend  to 
smooth  out  small  irregularities.  Compressive  forces  would  probably 
tend  to  decrease  the  movement  near  the  center  of  the  subsidence 
trough.  The  soil  probably  cracks  in  the  tension  zone  and  these  forces 
would  have  little  effect. 

The  movement  between  the  end  columns,  i.e.  the  lengthening 
or  shortening  of  the  overall  column  line  is  of  more  significance 
than  the  relative  movement  of  individual  columns.  Based  on 
measurements  taken  December  5,  which  are  believed  to  be  the  most 
accurate  due  to  the  care  in  correcting  for  the  temperature 
changes,  the  calculated  shortening  of  both  column  lines  G  and  H 
were  within  25  percent  of  the  measured  value.  Of  the  three 
separate  measurements  taken  on  each  column  line,  the  one  taken 
December  5  gave  the  largest  deviation  from  the  calculated  value. 

It  was  interesting  to  note  that  the  shortening  as  computed  from  the 
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other  two  sets  of  measurements  varies  by  10  to  15  percent  from  the  short¬ 
ening  calculated  on  the  basis  of  those  taken  December  5th.  More  than 
average  care  and  caution  was  used  when  taking  all  measurements  yet  the 
difference  between  actual  and  the  theoretical  shortening  is  not 
appreciably  different  from  the  difference  between  individual  measure¬ 
ments.  It  should  be  said  that  the  actual  measured  overall  distances 
agreed  within  0.05  feet  in  the  1097  feet  or  to  an  accuracy  of 
1  :  22,000. 

There  was  no  significant  difference  between  the  observed  or  com¬ 
puted  movements  of  column  lines  G  and  H  and  they  will  not  be  discussed 
separately.  Many  unknowns  existed  in  both  measured  and  computed 
horizontal  shortening  of  these  column  lines.  In  view  of  this  fact  it 
would  be  a  useful  contribution  to  find  an  explanation  for  the  pattern 
of  movement  alone.  The  close  agreement  of  the  actual  values  not  only 
supports  the  validity  of  the  basic  equation,  but  indicated  that  the 
results  obtained  from  its  use  are  of  the  correct  order  of  magnitude. 

Gravel  as  a  Beam 

The  preceeding  theoretical  computation  was  based  on  the  fact  that 
there  would  be  zero  lateral  movement  at  the  gravel  silt  interface. 

Another  similar  approach  was  to  consider  the  gravel  layer  as  a  beam  In 
deflection.  According  to  beam  theory  (29)  a  surface  exists  within  a 
deflecting  beam  which  suffers  no  lateral  expansion  or  contraction. 

This  is  called  the  neutral  surface  and  for  homogeneous  beams  of 
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symmetrical  cross  section  is  at  midheight.  For  small  deflections  the 
horizontal  movement,  v,  of  a  point  not  on  the  neutral  surface  relative 
to  another  point  the  same  distance  from  the  neutral  surface  can  be 
expressed  by  the  equation: 

V  s  y<*  (8) 

where  y  is  the  distance  from  the  neutral  surface  to  the  point  and  oc 
is  the  angle  change  at  the  neutral  surface. 

If  the  upper  gravel  layer  is  considered  as  a  deflecting  beam 
then  by  this  theory  y  becomes  equal  to  h/2  and  the  horizontal  displace¬ 
ment,  v,  of  a  point  at  the  surface  is  expressed  by: 

v  =  h  <*  (9) 

2 

This  is  exactly  one  half  the  value  of  the  previous  method,  or  equation  7* 
Adopting  this  formula  to  the  case  of  column  lines  G  and  H  would 
in  effect  reduce  the  computed  lateral  movement  at  all  points  by  one 
half.  As  shown  on  Figures  66  and  6?  this  results  in  close  agreement 
with  observed  data  in  the  central  zone,  but  the  disagreement  at  the 
ends  of  the  building  is  greater  than  when  computed  by  equation  7* 

Comparison  to  Movement  from  Stress  Computations 

According  to  the  beam  theory,  lateral  movement  would  occur  at 
both  the  upper  and  lower  surfaces  of  the  gravel  layer.  The  direction 
of  movement  at  the  upper  surface  would  be  toward  the  zone  of  maximum 
settlement  and  the  movement  of  the  lower  surface  would  be  equal  and  in 
the  opposite  direction  to  that  at  the  upper  surface.  Since  no  slipping 
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can  occur  at  the  gravel-silt  interface  the  beam  theory  assumption 
requires  lateral  outward  movement  of  the  silt  layer  also*  This  move¬ 
ment  at  the  lower  gravel-silt  interface  would  be  a  maximum  in  the  zone 
of  maximum  differential  settlement  and  decrease  to  zero  at  the  outer 
edges  and  at  the  point  of  maximum  settlement* 

In  one  respect  the  beam  theory  agreed  nicely  with  the  horizontal 
movements  computed  from  the  three  dimensional  stress  analysis  in 
Chapter  III*  Both  cases  showed  lateral  outward  movement  at  the  gravel- 
silt  interface*  This  outward  movement  was  of  approximately  the  same 
order  of  magnitude,  0.25  to  0*50  feet  maximum  for  both  cases*  An 
important  difference  is  the  fact  that  movement  based  on  beam  theory  was 
a  maximum  in  the  zone  of  maximum  differential  settlement  which  at 
Kitimat  was  near  the  middle  of  the  area.  The  movement  in  this  zone 
based  on  the  stress  analysis  was  a  minimum.  At  the  edges  where  the 
stress  consideration  gives  maximum  lateral  movement,  values  from  the 
beam  theory  are  small. 

Movement  computed  by  the  beam  theory  alone  gave  results  which 
agreed  well  with  the  observed  movement.  However,  if  movement  from  the 
beam  theory  were  modified  to  account  for  the  effect  of  a  non  vertical 
major  principal  stress  the  result  would  be  in  even  closer  agreement  to 
the  observed  condition.  Outward  movement  at  the  edges  based  on  the 
stress  consideration  would  tend  to  increase  the  total  outward  move¬ 
ment  at  the  gravel-silt  interface.  This  outward  movement  at  depth 
would  be  reflected  as  an  increased  inward  movement  at  the  surface  due 


..  •  *.  ••  o  i  •  o 

- 


r.  :■ 

* 

v  '  .  .  •;  '  ,  ■  'L 


* 

■  f  - 

'  .f  :■  ■'  ..  /,■ 

'  -  .  -  •  •  ( 

■  ■  ■  -  •  .  '  •  •  ;■ 

\  i  i  v  7.'  .  '  f  o:j  .r \:J  8 

o  .  )  '  -  : 


92 


to  rotation  of  vertical  elements  about  the  neutral  axis  of  the  gravel 
"beam”.  An  opposite  effect  would  occur  near  the  center  of  the  loaded 
area  where  lateral  movement  from  the  non-vertical  major  principal  stress 
is  small.  If  the  neutral  axis  were  considered  to  be  at  midheight  of  the 
gravel  layer  the  combined  movement  at  the  ends  of  the  potroom  would 
exceed  the  observed  movement  by  0.2  to  0*3  feet.  This  is  in  contrast  to 
the  fact  that  movement  due  to  beam  theory  alone  was  approximately  0.2  to 
0.3  feet  less  than  observed.  Although  it  was  not  shown,  it  is  probable 
that  some  of  the  movement  at  the  base  of  the  gravel  layer  due  to  a  non¬ 
vertical  major  principal  stress  would  be  damped  out.  This  would  mean 
that  less  movement  would  be  reflected  at  the  surface  and  the  net  result 
would  be  closer  agreement  with  the  observed  condition. 

It  seemed  a  significant  fact  that  according  to  the  available 
literature,  cases  of  horizontal  movement  associated  with  ground  sub¬ 
sidence  were  only  noticed  in  areas  where  the  seat  of  settlement  was  at 
some  considerable  depth  below  the  surface.  On  the  other  hand,  settle¬ 
ment  occurring  in  shallow  deposits  did  not  seem  to  be  bothered  with 
this  problem.  The  classical  example  of  ground  subsidence  is  at  Mexico 
Gity.  Here  the  entire  city  has  been  settling  as  a  result  of  lowering 
the  ground  water  table.  This  settlement  is  deep  seated  in  deposits  of 
volcanic  clay  which  extend  from  the  surface  to  an  unknown  depth. 

Cummings  (U)  and  Leonardo  (16)  have  discussed  the  foundation  problem  in 
this  city  with  respect  to  the  more  important  shallow  settlement  which 
results  from  heavy  building  loads.  Differential  settlements  from  this 
cause  have  amounted  to  as  much  as  $  or  6  feet  between  adjacent  buildings 
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yet  no  mention  is  made  of  any  horizontal  movement#  The  major  principal 
stress  within  the  soil  due  to  the  extreme  variation  in  building  loads 
would  almost  certainly  not  be  vertical.  The  fact  that  there  appears  to 
be  no  horizontal  displacements  resulting  from  this  shallow  seated  settle 
ment  helps  to  verify  the  statement  that  horizontal  movements  only  occur 
where  the  seat  of  settlement  is  at  some  distance  below-  the  ground  surfac 

General  Summary 

In  preceeding  sections  it  was  shown  that  the  distance  between 
the  ends  of  Potroom  2A  have  shortened  approximately  one  half  foot  since 
the  building  was  constructed,  it  was  further  shown  that  this  shortening 
has  occurred  mainly  in  the  region  of  maximum  differential  settlement. 
Theoretical  time  curve  matching  showed  that  no  consolidation  has  taken 
place  within  the  upper  gravel  layer  alter  completion  of  building  const¬ 
ruction,  hence  all  building  movement  has  been  a  result  of  consolidation 
within  the  lower  silts  and  clays. 

Subsidence  as  a  result  of  mining,  or  removal  of  oil  or  water 
from  deep  below  the  surface  has  resulted  in  horizontal  movements  of  the 
ground  surface  over  a  wide  area.  No  such  horizontal  movements  have 
been  recorded  for  settlement  due  to  shallow  seated  consolidation. 
Research  into  the  problem  of  mining  subsidence  has  established  the 
general  trend  of  this  lateral  displacement  to  be  toward  the  zone  of 
maximum  subsidence  with  the  greatest  movement  occurring  in  the  zone  <f 
maximum  differential  settlement.  Nhen  plotted  in  the  same  manner,  the 
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lateral  movement  along  Potroom  2A  agreed  in  general  with  that  of  mining 
subsidence* 

Two  equations  were  proposed  to  account  for  the  observed  horizon¬ 
tal  movement*  Both  equations  assumed  rotation  of  a  vertical  element 
about  a  neutral  axis  in  the  upper  gravel  layer#  They  were  in  fact  modi¬ 
fications  of  the  beam  theory  applied  to  the  upper  ground  crust*  When 
the  neutral  axis  was  assumed  to  be  at  the  base  of  the  gravel  stratum 
the  computed  horizontal  movements  were  too  large  in  the  region  of  maxi¬ 
mum  settlement,  and  slightly  too  small  toward  the  ends  of  the  potroom* 
When  the  neutral  axis  was  assumed  at  midheigho  of  the  gravel  layer  the 
computed  horizontal  movements  agreed  with  the  observations  In  the 
region  of  maximum  differential  settlement,  but  were  too  small  at  the 
ends  of  the  potroom. 

The  results  of  a  model  test  verified  the  findings,  both  computed 
and  observed,  that  there  was  movement  of  the  ends  of  the  column  lines 
toward  the  zone  of  maximum  settlement. 

An  analysis  was  also  made  on  the  effect  of  consolidation  proceed¬ 
ing  in  the  direction  of  the  major  principal  stress.  When  the  computed 
movements  based  on  beam  analogy  were  modified  to  account  for  the  effect 
of  a  non-vertical  major  principal  stress  the  results  showed  even  better 
agreement  with  the  observed  movement* 

It  should  be  noted  that  the  results  from  beam  analogy  computations 
alone  gave  results  which  were  in  good  agreement  with  actual  conditions. 
The  effect  of  modifying  the  computations  by  accounting  for  the  non- 
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vertical  major  principal  stress,  while  beneficial,  was  not  of  sufficient 
importance  for  it  to  be  a  routine  procedure  for  future  computations  of 
this  kind# 

Conclusions 

In  the  foregoing  discussion  it  was  shown  that  horizontal  movement 
along  the  column  lines  could  be  explained  on  the  basis  of  the  'Upper 
gravel  layer  acting  as  a  deflecting  beam  resting  on  deep  compressible 
material#  It  was  further  shown  that  the  effect  of  the  non-vertical 
major  principal  stress  within  the  compressible  material  acted  in  such 
a  way  as  to  cause  shortening  of  the  column  line,  then  the  results  of 
the  beam  analogy  were  modified  to  account  for  the  effect  of  a  non¬ 
vertical  major  principal  stress  the  resulting  computed  movement  was 
in  closer  agreement  to  the  observed  movement  than  that  computed  by  the 
beam  analogy  alone# 

Future  Movement 

In  Chapter  III  it  was  shown  that  the  settlement  within  the  silt 
and  clay  layers  was  at  the  time  of  writing  approximately  60  percent 
completed.  Since  the  column  lines  have  shortened  one  half  foot  at 
the  time  of  writing  it  is  probable  that  the  future  shortening  will 
amount  to  approximately  0#3  feet  and  the  total  ultimate  shortening  will 
be  approximately  0.8  feet  when  the  settlement  is  completed. 
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Damage  to  the  Supers  true ture 

A  general  observation  of  the  distortion  of  the  superstructure 
indicates  the  source  of  damage  to  be  the  lateral  movement  at  the 
ground  surface.  A  detailed  analysis  of  these  distortions  is  beyond 
the  scope  of  this  study.  Should  research  along  this  line  be  attempted 
in  the  future  the  pertinent  data  can  be  found  on  Figures  68  and  69 
enclosed  herewith. 
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APPENDIX  A 

TIME  FACTOR  VALUES 


The  following  are  values  for  time  factor  (t)  at  different  degrees 
of  consolidation  (U),  which  were  used  in  this  study  to  match  a  theoreti¬ 
cal  time  curve  to  the  observed  settlement  time  plots  for  Settlement 
Gauges  1  to  15*  These  values  were  taken  from  Curve  I  of  Reference 
2k}  page  237  and  represent  the  case  of  linear  variation  of  initial 
excess  hydrostatic  pressure* 


u$ 

T 

10 

.008 

20 

.031 

30 

.071 

ho 

.126 

5o 

.197 

60 

.287 

70 

.U03 

80 

.567 

50 

.8it8 

Where  required,  time  factors  for  intermediate  values  of  U  were 
computed  from  the  approximate  expressions 
U  less  than  60$  T  =  7T 

~"IT" 


U  greater  than  60$  T  =  -0*9332  log^g  (1-U)  -0.0851 
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Symbol  Name  Dimension 

z ,x,y  Rectangular  or  Cartesian  coordinates  ft. 

z  positive  down 

x  positive  south 

y  positive  west 

°C  Angle  of  inclination  between  the  subsidence 

curve  and  original  ground  surface 
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Sign  Convention 
AN  jo  NOTATION  for 

Axes  and  Stresses 


OC  -  Vertical  Compo/renf 
L/J  03  / OD  -  Hot^iz  .  Corr/ps. 


o  /  OP3  Coi  <*  8 


08  -0/4  Cos  cX. 

03  -  a%o$ 

•  oS  -  —  *  Cosox 
^OS'X 


S/rti/ori/y 

OD-  OC 


Cosfi 

£*sy 


O/r'e c.i~i o n  of  PnwapoJ 
S+ress  CT 


Three  Dimensional  Vectors 


Solve  the  following  formula  f or  principal  stress  -  G 

c r3-[c;  +  <r^ci)cz  +  fe  +  cr^  +cryci-tc;cx  <r 

'  \p%  (7^  Gz"  +  2.  Zy£  7xz  Ttfy  TT  2yz  —  -  y  -  *z  “  CT  Xyz  }  -  0 


Solve  -for  the  direction  cosines 

2 


ffCrrr- 

Cos'g  +  Cos4  y  =  1 

and  an^  two  of  the  -following- 


Co  sZp< 


(cr  -G~ 

)  Cos 

+  Txy 

Gos  (p>  A-  Ticz  C*S  ^ 

=  0 

ry* 

^05  c< 

-  (q-<r) 

CoS  (3  -r  T-yH  Co%Y 

=  0 

r2* 

Cos  oC 

+  'Tyz 

Cos  §  +  (<^-r)  Cos  y 

=0 

5ol\/e  -for  -the  horizontal  component-  in  terms  of 
the  yer-pieal  Componenf:  _ 

Horn- Comp  (x.d^ecuan)  =?  Coi^Z  *  Vert.  Comp. 

Horn.  Comp /v-a.v«h»ni  =  *  Ver-h  Comp. 
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HORIZONTAL  DISPLACEMENT  CAUSED  BY  MINE  SUBSIDENCE 


Basic  Assumptions  Displacement  Equations 


v  =  horizontal  displacement  at  the  surface  of  any  point  P  (  always  toward  center) 
x  =  horizontal  distance  from  center  of  trough  to  point  P 
y  =  subsidence  of  point  P 

m  =  total  maximum  subsidence  at  the  center  of  the  trough  =  thickness  of  seamm 
F  =  area  between  original  ground  and  subsidence  curve  from  center  to  point  P 
R  =  distance  from  center  of  subsidence  trough  to  edge  of  excavation  (point  of 
maximum  slope  of  subsidence  curve) 
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SUMMARY  OF  OBSERVATIONS  OF  MINE  SUBSIDENCE  CURVES 

1#  The  maximum  sinking  is  in  the  middle  of  the  trough* 
Here  there  is  no  shifting* 

2*  All  the  shifting  is  directed  to  the  center  of  the 
excavation.  The  maximum  shift  is  at  that  point 
where  no  tension  and  compression  occurs;  here  is 
the  point  without  stress  and  the  changing  point 
of  curvature  as  well  as  the  steepest  dipping  in 
the  course  of  the  trough. 

3*  The  region  with  the  maximum  tension  lies  near  the 
outer  edge  of  the  trough* 

Ij..  The  region  with  the  maximum  compression  lies  near 
the  inner  edges  of  the  trough* 

The  sum  of  the  tensions  is  equal  to  the  sum  of 
the  compressions,  the  total  sum  being  zero. 


(After  Rellensmann) 
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